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Abstract 

Presented is description of kinematics and dynamics of material points 
with internal degrees of freedom moving in a Riemannian manifold. The 
models of internal degrees of freedom we concentrate on are based on the 
orthogonal and affine groups. Roughly speaking, we consider infinitesimal 
gyroscopes and homogeneously deformable gyroscopes (affienly-rigid bod- 
ies) in curved manifolds. We follow our earlier models of extended rigid 
and afRnely-rigid bodies moving in a flat space. It is well known that 
in curved spaces in general there is no well-defined concept of extended 
rigid or afhnely-rigid body. Our infinitesimal models are mathematically 
well defined and physically they may be interpreted as an approximate 
description of "small" rigid and afRnely-rigid bodies. We derive equations 
of motion and show how internal degrees of freedom interact with spa- 
tial geometry, first of all with the curvature but also with the torsion. 
Integrability and degeneracy problems are discussed. 

Keywords: internal degrees of freedom, Riemannian manifolds, infinitesimal 
gyroscopes and affinely-rigid bodies, coupling of space geometry and internal 
degrees of freedom. 



1 Introduction 

The primary concept of Newton mechanics is that of the material point mov- 
ing in three-dimensional Euclidean space. A good deal of the theory depends 
only on the afhne sector of geometry. The metric structure becomes essential 
when constructing particular functional models of forces; the concepts of energy, 
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work, and power (time rate of work) also depend in an essential way on the met- 
ric tensor. The Galilei relativity principle implies that, as a matter of fact, it 
is not three-dimensional Euclidean space but rather four- dimensional Galilean 
space-time that is a proper arena of mechanics. This space-time has relatively 
complicated structure, does not carry any natural four-dimensional metric ten- 
sor and fails to be the Gartesian product of space and time. There exists the 
absolute time, but the absolute space docs not. In the sequel we concentrate on 
the other kind of problems, so the analysis of the subtle space-time aspects will 
be almost absent in our treatment. Newton theory becomes essentially realistic 
and viable when multiparticle systems arc analyzed. It is just there where met- 
rical concepts become almost unavoidable, because it is practically impossible 
to construct any realistic model of interparticle forces without the explicit use 
of the metric tensor. Extended bodies are described as discrete or continuous 
systems of material points. Their motion consists of that of the center of mass, 
i.e., translational motion and the relative motion of constituents with respect 
to the center of mass. The total configuration space may be identified with 
the Cartesian product of the physical space (translational motion) and the con- 
figuration space of relative motion. In many physical problems the structure 
of mutual interactions leads to certain hierarchy of degrees of freedom of the 
relative motion; in particular, some constraints may appear. The effective con- 
figuration space becomes then the Gartesian product of the physical space and 
some manifold of additional degrees of freedom. There arc situations when this 
auxiliary manifold and the corresponding dynamics are postulated as something 
rather primary then derived from the multiparticle models. Usually the guiding 
hints are based on some symmetry principles. In this way the concept of internal 
degrees of freedom replaces that of relative motion. Sometimes it is a merely 
convenient procedure, but one can also admit something like essentially internal 
degrees of freedom not derivable from any multiparticle model. After all, the 
very concept of the material point is an abstraction of a small piece of matter. 
Why to reject a priori an abstraction of a small and particularly shaped piece 
of matter, thus, the material point with extra attached geometric objects, as 
something primary for mechanics? In this way the configuration space becomes 
as follows: 

Q = M X Q"*, 

the Cartesian product of physical space M (translational motion) and some 
internal configuration space Qint- The usual d'Alembert procedure of deriving 
equations of motion is not then reliable, perhaps just essentially inadequate and 
should be replaced by something else, probably based on appropriate invariance 
assumptions. 

In standard quantum theory there are quantities which have the status of 
essentially internal variables, e.g., spin. If M is not the Euclidean space but some 
general manifold, then the concept of essentially internal degrees of freedom 
becomes even more justified, as we shall see. 

And just in connection with this, the last step of generalization: why Q = 
M X <5'"*? Perhaps any spatial point a; G M has its own manifold of internal 
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variables Q™*? Then the total configuration space would be 

Q=[j 

xeM 

III this way, the conc;ept of fibre biuidle appears as a most natural mathematical 
framework for describing internal degrees of freedom. 

As mentioned, such concepts are particularly natural when M is a general 
manifold endowed with some geometric structure based on the afhne connec- 
tion, metric tensor, or both (interrelated or not). In any case, the primary 
mathematical concept underlying any model of internal degrees of freedom is a 
principal fibre bimdle (Q.M, tt), where M denotes the base manifold (physical 
space), Q is the total bundle manifold (configuration space), and n : Q ^ M is 
the bundle projection. We shall often use the following abbreviation for fibres: 
Qx = iT~^{x). Obviously, in the relativistie theory a more adequate formulation 
is one using the space-time manifold as a basis of the fibre bundle. The same is 
true in non-relativistic theory when the problems of Galilean relativity principle 
are important. However, in this treatise we do not touch such problems or do 
it merely exceptionally. 

As is well known, the general fibre bundle (Q,M, tt) does not need to be 
diffeomorphic with the Cartesian product M x Q'"*. And even if such a dif- 
feomorphism does exist it is not in general unique in the sense that there is 
no canonically distinguished choice (there are exceptions like, e.g., tangent and 
cotangent bundles over Lie groups, bundles of linear frames or co-frames over 
Lie groups, and so on). Therefore, in general, the total motion g -.R^ Q does 
not split in a well-defined way into translational and internal motion. More 
precisely, only translational motion is well defined as a projection 

g^^ := w o Q -.R^ M. 

Similarly, generalized velocity g{t) G Tg(t)Q does not split into translational and 
internal velocities. Only the first one is well defined as a T7r-projection of g{t) 
to the tangent space T^(g(t))M; obviously, it is identical with Quit) e Tg^^^i^M. 
Without additional geometric objects the time-rate of internal configuration is 
not well defined. 
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Fig. 1 

The very concept of internal motion and internal trajectory is meaningful 
only when Q = M x Q'"*, i.e., when Q is the Cartesian product. In general, 
when the bundle {Q,M,tt) is nontrivial (perhaps even not trivializable) and no 
additional connection-like structures in Q are defined, it is only the total motion 
in Q and the translational one in M that are well defined. The same concerns 
velocities. 

In original Newton theory M is the Euclidean three-dimensional space. We 
shall consider the manifolds with a more general geometric structure and the 
fibre bundles (Q, M, tt) over them. Having in view both the mathematical clarity 
and certain nonstandard physical applications we admit the general dimension, 
i.e., dim M = n. 

2 Kinematics and canonical formalism for affine 
models 

In a long series of earlier papers we have developed the theory of extended 
affinely-rigid bodies in Euclidean spaces [Tl[Il[Til[Ill[ni[Il[Ml[151lli[ll 
[2S[3ll[3llM|3l[3i|3Zl|3i|4ll|43|42|4i|43 Strictly 
speaking, kinematics of such objects preassunies only affine structure of the 
physical space. Affine constraints mean that all affine relations between mate- 
rial points remain invariant during admissible motions; material straight lines 
continue to be straight lines, their parallelism is preserved, etc. Assuming that 
all metrical relations are preserved one obtains the rigid body in the usual sense. 

The concepts of extended metrically- and affinely-rigid bodies break down 
when the Euclidean or affine space is replaced by a differential manifold with 
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geometry given by the metric tensor, afRne connection, or both of them (in- 
terrelated or not). For example, let {M,g) be a Riemann space, where M is a 
manifold and is a metric tensor defined on it. An extended continuous system 
of material points moving in M is metrically rigid if all infinitesimal distances 
are invariant during any admissible motion. Of course, then also finite along- 
geodetic distances are constant. If some standard reference configuration is 
fixed, the configuration space becomes identified with Diff(M, g), the isomctry 
group, i.e., the group of all diffeomorphisms (p : M ^ M preserving the metric 
tensor, (f*g = g. But it is well known that for the generic Riemann space (M, g) 
with not vanishing curvature tensor R[g\ 7^ it is rather typical that the isom- 
etry group is trivial, i.e., Diff(M, (7) = {idM}- Obviously, the highest possible 
dimension of Diff(M, g) in an n-dimensional M equals n(n + l)/2. This highest 
dimension is attained only in constant-curvature spaces, like, e.g., spheres and 
pseudospheres in M"+^. And in fact, in such manifolds the concept of a finite 
extended rigid body continues to be well defined. But such manifolds, with 
n{n + l)/2 independent Killing vectors, are extremely exceptional in category 
of all Riemann manifolds. Let us also mention that even if R[g\ = 0, i.e., (M, g) 
is locally Euclidean, but M is topologically not equivalent to R" (e.g., torus 
T" = R^/Z" with the metric tensor inherited from M"), there may be obstacles 
against the existence of the global n{n+ l)/2-dimensional isometry groups. So, 
in a generic Riemann manifold the concept of extended rigid body fails to be 
well defined. 

The same concerns extended affinely-rigid bodies. Let (M, T) be an affine- 
connection space, i.e., a differential manifold M endowed with an affine con- 
nection r. This connection may be quite arbitrary, not necessarily one derived 
explicitly from some metric tensor g (at this stage the metric structure does not 
need to exist at all). Let V denote the covariant differentiation corresponding 
to r. We say that € Diff(M) is an afhne transformation of (Af, F), i.e., that 
G Diff(M, F) if for any pair of vector fields X,Y on M the following holds: 

V^.x {^*Y) = ip^ {VxY) , 

where, as usual, V x denotes the covariant differentiation along the vector field 
X. Roughly speaking, the above property means that the V-operation is "trans- 
parent" with respect to the transformation ip. In the usual flat affine space this 
is just an equivalent definition of the n(n+ l)-dimensional affine group GAf(A'/). 
However, when the curvature tensor 7?,[F] of the affine connection F is not van- 
ishing, then, as a rule, the dimension of Diff(M, F) is smaller than n{n + 1) and 
the generic situation is that Diff(Af, F) = {idAf}, i.e., the only affine transfor- 
mation is the identity mapping. Therefore, in principle the concept of extended 
affinely-rigid body breaks down. 

However, "very small" regions of (M, F) approximately look like flat affine 
spaces; similarly, "very small" regions of Riemannian manifolds (M, g) look 
like Euclidean spaces. Therefore, approximately, one can consider "very small" 
rigid and affinely-rigid bodies in non-Euclidean manifolds. This is always an 
approximation, the better one, the more body shrinks. Everything becomes 
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rigorous in the limit of vanishing size. Roughly speaking, the body is not any 
longer injected in the physical space M, but in a tangent space T^M, where 
X G M represents the spatial position of the body "as a whole" and is a remnant 
of the centre of mass position in the flat-space theory. And configurations of 
affine bodies in the linear space T^M may be identified with linear frames, i.e., 
ordered bases in T^M. 

This is a natural analogue of some description used in mechanics of extended 
affine bodies. Namely, let V be the linear space of translations in the physical 
affine space M. The configuration space of extended aflanely-rigid body in M 
may be identified with M x F{V), where F{V) denotes the manifold of linear 
frames in V. This is meant in the following sense: a; € M is an instantaneous 
position of the centre of mass, and the frame e = (ei, . . . , e^, . . . , e„) e F{V) 
is materially frozen into the body, i.e., co-moving with it. More precisely, if 
a = (a^, . . . , a") are Lagrange (reference) coordinates of some material point (its 
identification labels), and configuration is given hy q = {x, e) = {x;ei, . . . , e„) e 
Q = M X F{V), then the current spatial position y e M of this a-th point 
satisfies 

xy = a^BK, 

where xy denotes the radius- vector of y with respect to the instantaneous posi- 
tion a; e M of the centre of mass. Analytically, 

y' =x' + e^Ka^ 

with respect to some Cartesian coordinates. 

Essentially the same remains true for infinitesimal bodies injected in tangent 
spaces TxM. They are somewhere placed in space, x G M, and have the extra 
attached internal variables ex € T^M, K = l,n. Let us describe it in geometric 
terms. 

The configuration space of infinitesimal affinely-rigid body moving in the 
physical space M is given by the manifold FM of linear frames in M, 

Q = FM= U F,M, 
xeM 

where F^M denotes the manifold of linear frames in the tangent space T^M. 
Obviously, just as F{V) is an open submanifold of 

y» = y X F X • • • X y, 
^ ^ 

n terms 

SO FM is an open submanifold in the Whitney sum of fibre bundles 
0TM= IJ T^MxT^M X ■■■ x T^M^ . 

These open subsets consist of linearly independent linear n-tuples. 

To be more precise, in the theory of the flat-space extended continuous bodies 
one should not use the total F{V), but rather one of its connected components. 
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e.g., one positively oriented with respect to some fixed standard orientation. 
And so is in mechanics of structured material points; by the way, it is tacitly as- 
sumed that M is orientable. Although, one can argue that if once the framework 
of classically imaginable extended bodies is left and affine degrees of freedom 
are essentially internal, the exotic orientation-changing motions might be per- 
haps admissible and even (at least mathematically) interesting. The problem 
becomes particularly interesting, perhaps just mathematically exciting, when 
M is not orientable, like, e.g., Mobius band. 

The manifold FM carries a natural structure of the principal fibre bundle 
over the base M I53j . The projection tt : FM M assigns to any linear 
frame the point at which it is attached, thus tt [F^M) — x. The structural 
group GL(n,R) acts on FM according to the standard rule; thus, for any L e 
GL(n, R) we have that 

FM3e={..., eA, ...) eL :=(... , e^L^A, • ■ •). (1) 

Obviously, dim FAI — n{n + \); this is the number of degrees of freedom (n 
translational and internal ones). For any linear frame e there exists a unique 
dual co-frame e — (. . . , e^, . . .), where 

and for any x € M, v € T^M and p S T*M the symbol (p, v) denotes the 
evaluation of the linear function p on the vector v. 
The manifold of all co- frames, i.e., 

Q* ^ F*M = IJ f;m, 

xeM 

is canonically diffeomorphic with FM just in the sense of duality. Therefore, 
the configuration space of infinitesimal affine body may be represented either 
as FM or F*M; it depends on the particular problem which representation is 
more convenient. Just as previously, F*M is an open subset of the Whitney 
sum 

T*M ^ IJ T*M X T*M X • • • X T*M; 

n terms 

it consists of linearly independent rt-tuples. The natural projection onto the 
base manifold will be denoted by tt* : F*M — + M, and then tt* {F*M) — x. 
Obviously, F*M also carries the structure of the principal fibre bundle with 
GL(n, R) as a structural group: 

F*M3e^ (...,e^,...) ^eL:= (. . . , L-^-^^e^, . . .) (2) 

for any L e GL(n,R). 

There is one subtle point, namely, GL(7i, R) is not connected. It con- 
sists of two connected components, i.e., the subgroup GL+(n,M) of positive- 
determinant matrices and GL~ (n, R) is the coset (not a subgroup, of course) 
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of negative-determinant matrices. In theory of Lie groups, fibre bundles and 
connections one usually deals with connected groups. Thus, it seems rather 
natural to deal with reductions to the subgroup GL+(n,M). But this has again 
to do with orientability of M. Then the unconnected manifold FM is replaced 
by M , i.e., the manifold of linear frames positively oriented with respect to 
some fixed orientation on M. 

Any system of local coordinates x*, i = 1, n, on M induces local coordinates 
(x*, e*A), A = l,n, on FM and (x% e^i) on F*M, where e*^, e^i are respectively 
components of ca, with respect to coordinates x% and 

To avoid the crowd of symbols we do not distinguish graphically between and 
their pull-backs to FM and F*M. 

In certain problems it is convenient to admit the local, x-dependent action 
of the structural group, like in gauge theories. This means that we consider 
(sufficiently smooth) fields L : M — > GL(?t., R); they form an infinite-dimensional 
group under the pointwise multiplication, 

{LiL2){x) = Li{x)L2{x). 

And this group acts, also in a pointwise way, on FM and F*M. Thus, if 
e e FxM, e e F*M, then the action of L is given by 

ei~*eL{x), ei~+eL{x). (3) 

From the mechanical point of view the structural action of GL(n, M) on Q = 
FM and Q* — F*M corresponds to material transformations. We shall use 
the term "micromaterial transformations" . This is exactly the infinitesimal 
limit of the usual material transformations. In fact, material points with affine 
internal degrees of freedom may be interpreted as affinely-rigid bodies injected 
into instantaneous tangent spaces T^M, where x G M is the spatial position of 
the body. But it is obvious that e G F^M is canonically identical with some 
linear isomorphism of M" onto T^M (similarly, e S F^M is a linear isomorphism 
of TxM onto K"). In this way, R" plays the role of the "micromaterial space" 
(Lagrange variables) and T^M is the "microphysical space" (Euler variables) of 
infinitesimal affinely-rigid body. The frame e € F^M, i.e., the "placement" is 
then exactly the counterpart of what was denoted hy (p G LI(?7, V) in our papers 
about extended fiat-space bodies [3T1I31IS3I31IM1IS71I31I1I1I11I11I151I1S1 
HzlllHllli [501151]. 

And what does correspond to spatial or physical transformations, i.e., what 
are "microspatial" kinematical symmetries? Here the situation is more compli- 
cated. Obviously, microspatial transformations must act on the left in tangent 
spaces TxM . But for different points x G M the tangent spaces TxM are logi- 
cally different linear spaces and their linear groups GL(r2;M) are also logically 
different sets. Without additional strong geometric structures there is no natu- 
ral isomorphism between different GIj{TxM). So, there is only one possibility. 
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in an essential way infinite-dimensional. Namely, the fields T of mixed non- 
degenerate second-order tensors on M ^ 

MBx^T^e GL(r,Af) C L(r,M) ~ T\{T^M), 

give rise to the following transformations of FM , F*M: 

F^M 9 e = (..., CA, .. .) ^ (. . . , o eA, . . .) , (4) 
F:M3e={...,e^,...) ^ {. . . , o T-\ . . .) . (5) 

In a structure-less manifold M and in a connection manifold (M, F) with not 
vanishing curvature tensor 7?,(r), there are no natural isomorphisms between 
different T^M, therefore, nothing like the "constancy" of T may be defined. If 
the components field T'j is accidentally constant in some coordinates a;*, then 
in other coordinates it is no longer true. Therefore, by its very nature the above 
group of left-acting transformations is infinite-dimensional, fiexible. Of course, 
in the right-hand side action of GL(n,M) ^ L may be put x-dependent like 
in gauge theories of fields and continua, but need not to be so; its constancy is 
well defined. L € GL(n,]R) are just matrices by their very nature, not matrix 
representants of linear mappings in T^M with respect to some coordinates. 

When some coordinates x'' are fixed in M, the second-order mixed tensor 
field T is represented by a system of functions, i.e., components T*j(a;°), 

T = T',{x)^(g>dx^. 

Analytically, the action of T on the configuration space Q is described by 

(...,x^...;...,e'A,...) ^ (...,a:^...;...,r,(x)e^■A,...). (6) 

Similarly, x-dependent L e GL(n,R) act as follows: 

(...,x^...;...,e^A,...) ^ (...,x^...;...,e^si^^A(x),...). (7) 

In particular, this formula holds for global, a;-independent L E GL(n,R). It 
is seen that when x G M is kept fixed, these are just the well-known formulas 
for an affinely-rigid body in a fiat space T^M. The difference is that there is 
no counterpart of translations and general affine mappings in M and in the 
total FM. The reason is that, as mentioned, even if M is endowed with some 
affine connection F with not vanishing curvature tensor, the group of affine 
transformations is generically trivial or its dimension is smaller than n{n -1-1). 
And if M is completely amorphous, i.e., even if any affine connection is not 
fixed, then it is only the total diffeomorphism group Diff(M) that may replace 
the group of spatial affine transformations in mechanics of extended affine bodies 
in a fiat space. 

In mechanics of affine extended bodies [STlISlISaiMlISaEZlIMlllIlSlIll 
l45l l46l l47l l48l [49l [50l [5T] we considered various kinds of additional constraints; 
for obvious reasons the most important of them were metrical constraints, i.e., 
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rigid body in the literal sense. Not only affine relations between material points 
are then preserved but also metrical ones, i.e., distances and angles. In flat-space 
theory one is dealing then, also on the kinematical level, with the Euclidean 
structure (M, V,g), where V is the linear space of translations in M, and g £ 
V* ^ V* is the metric tensor of M. The configuration space of the rigid body 
may be identified with Q ~ AI x F{V, g), where F{V, g) C Fiy) is the manifold 
of g-orthonormal frames e = (ei, . . . , ca, . • ■ , e„), i.e., 

gieA,eB) = gijC^AC^ B = Sab- 

More precisely, one should use rather some connected component of F{V,g), 
e.g., the manifold F~^(V,g) of g-orthonormal frames positively oriented with 
respect to some fixed orientation of M. 

When (M, g) is a general Riemann manifold, then, as a rule, there are no 
extended rigid bodies with the usual number of n{n + l)/2 degrees of freedom; 
the only exception are constant-curvature spaces. In a generic case the not 
trivial isometrics do not exist at all. But just as in afhne theory we can speak 
about infinitesimal rigid bodies. They describe in a good approximation the 
behaviour of very small "almost rigid bodies" , in which the changes of distances 
between neighbouring particles are higher-order small in comparison with the 
initial distances themselves. 

The configuration space of infinitesimal rigid body in (Af , g) may be identi- 
fied with F{M,g), i.e., the manifold of all g-orthonormal frames in all tangent 
spaces of M. Again, to be more precise, if M is orientable, we should restrict 
ourselves to F'^{M, g), i.e., the connected manifold of g-orthonormal frames pos- 
itively oriented with respect to some fixed orientation in M. Obviously, F{M, g), 
F~^{M,g) are n(n -I- l)/2-dimensional manifolds; there are n translational de- 
grees of freedom (Af) and n(n — l)/2 rotational ones (fibres {FxM, gx))- Just 
as in affine model, it does not matter whether we use F{M, g) or the manifold 
F* {M, g) of g-orthonormal co-frames. Projections onto the base manifold M are 
just the restrictions of the previous tt, tt* to submanifolds F{M,g), F*{M,g); 
for brevity we denote them by the same symbols. The structure groups of 
F{M,g), F+[M,g) are respectively 0(n,R) C GL(n,M), SO(n,M) C GL(n,R); 
they act on the bundle manifolds on the right, just in the sense of formulas ([T]), 
Similarly, the left-hand-side spatial transformations are given by (jH) and 
([5]), where now the field T takes values in orthogonal groups of {TxM,gx), i.e., 

M3X ^ e 0(T,A/,g,) C GL(r,A/), 
M3X ^ TxeSO{TxM,gx)cO{TxM,gx); 

obviously, gx E T*M ®T*M is the metric tensor oiTxM , i.e., the value of the 
field 5 at a; e Af . 

What concerns the very description of degrees of freedom, infinitesimal 
affinely-rigid body is well defined in any quite amorphous differential mani- 
fold. Obviously, the usual, i.e., metrically- rigid body is meaningful only when 
M is endowed with some metric tensor g. 
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But, as mentioned, even on the purely kinematical level there are some 
difficulties with systems with internal degrees of freedom. Namely, for any 
motion g : R ^ FM it is only the tangent vector g{t) G Tg(^t)FM and its 
projection £>tr(i) € T^(^g(^t^-)M (see Fig. 2) that are well defined velocities, i.e., 
respectively the total generalized and translational velocities. 

There is no well-defined time rate of internal variables evolution. The mini- 
mal geometric structure necessary to define it is a connection on the fibre bundle 
{Q,M,tt). In mechanics of infinitesimal affine bodies and infinitesimal gyro- 
scopes this is simply the affine connection. Therefore, from now on we assume 
that some affine (more precisely linear) connection F is fixed in M. Analytically 
it is given by the system of components T^jk transforming under the change of 
coordinates on M according to the well-known linear-inhomogeneous rule. 
In modern differential geometry [19[ 153] linear connection is described by some 
differential one-form u on FM (or F{M,g)) taking values in the Lie algebra of 
the structural group GL(n,]R)' ~ L(n,R) (or SO(n,R)', i.e., the space of skew- 
symmetric matrices). It satisfies certain conditions that are quoted in [191 153] . 
Analytically uj is represented by the system of differential one-forms cu^ l related 
in the following way to coordinate-dependent quantities T^jk- 

OJ^B = e^^ {de^B + T'jke^Bdx'') . 




One uses also the connection-independent canonical R" -valued differential 
one-form 9 on FM; it is intrinsically defined and its coordinate representation 
is given by 
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It vanishes when evaluated on vectors tangent to the fibres F^M, Fx{M, g), i.e., 
vertical ones. The kernel of oje is transversal to the space of vertical vectors; 
its elements are by definition "horizontal" vectors in TeFM, TeF{M, g). The 
subspaces of vertical and horizontal vectors at e are denoted respectively by Ve , 
He- Obviously, T^FM = Ve®He, and similarly for F{M,g). Roughly speaking, 
horizontal vectors establish isomorphisms between fibres over infinitesimally re- 
mote points of M. Therefore, infinitesimally, the concepts like the change of 
internal state and the velocity of internal motion become meaningful. 

At any e, the co- vectors LUe,de & T*FM {T*F{M, g) in the gyroscopic case) 
form a basis of T*FM. The dual basis in T^FM {TeF{M, g)) consists of vectors 
denoted by E^b, FIa and 

{uj'^l.E^b) = <5^s<5^L, (lo'^l^Ha) - 0, 

=0, =<5^A. 

One can easily show that, after identifying vector fields with differential opera- 
tors [ini [53] , we have that 

E L=eL— — , iii = e L — - r j^e-'A-^-T^ • 

At any point e, the corresponding vectors [E^ ^ are vertical, i.e., tangent to 
the fibres, {E^ l) ^ e Ve, whereas {Hk)^ are horizontal, {Hk)^ G He- Dually to 
the situation with co-vectors w^Sj ■, now E^ are connection-independent, 
whereas depend explicitly on the connection. In the literature [121 [S3] Hk 
are referred to as standard horizontal vector fields and E^ l as fundamental 
vector fields. E^ l are infinitesimal generators of the action of the structural 
group GL(n,R) on FM. Roughly speaking, Hk generate the parallel transport 
in the sense of F. 

The splitting T^FM = Ve ® Hg enables one to decompose every vector at e 
into vertical and horizontal parts. In particular, this may be done for generalized 
velocities g{t) e Tg(j.-^FM . Projecting the horizontal component to M, we do 
not obtain anything new, just the translational velocity Qttit) G T^{g{t))M ■ The 
vertical part, however, is a new, F-depending object. It is a measure of the 
internal velocity, a kind of the time rate of internal configuration. And, as 
expected, it simply coincides with the system of covariant derivatives of the 
frame vectors ca along the curve describing translational motion. If x*, a are 
coordinates on FM induced by those on M, and if instead of sophisticated 
symbols o p, c'a o g we simply write x'^{t), e*A(i) for the time dependence 
of generalized coordinates of the object, then translational velocity v has the 
components 

^ dx' 

whereas the internal velocity is given as follows: 
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It is very important to stress that if (M, F) is non-Euchdean, i.e., the curvature 
and torsion tensors do not vanish: 

T^^jki = dkT"^ ji — diVji- + r\i.^°'ji — r\iT°'ji~ ^ 0, 

then the system (u% V^a) is the aholonomic velocity in the sense that there are 
no coordinates g', q'^A in FM for which the following could hold: 

dq^ ■ dq'A 

Generalized velocities at e € ty~^{x) arc represented by the (n + l)-tuplcs of 
vectors {V\ . . . ,Va, ■ ■ ■) & (T^M)"^^ attached at a; = 7r(e) G M. Therefore, the 
afHne connection T enables one to identify the tangent bundle TFM with the 
following fibre bundle over M: 



•,2n+l 



It is an open subset of the last Whitney sum; namely, the first n-tuplc is the 
submanifold of {TxM)" consisting of linearly independent systems. Obviously, 
for any fixed afiine connection T the above-mentioned diffeomorphism of TFM 
onto TFM is canonicial. 

Similarly, canonical momenta at e G 7r~^(a;) are represented by the (n + 1)- 
tuples of co-vectors at x = 7r(e), i.e., (P; . . . , P^, . . .) e (T*M)"'+\ 

It is seen again that F establishes a distinguished diffeomorphism of the 
cotangent bundle T*FM (phase space of the system) with the bundle 



T*FM = y F^M e (t;m)"+^ c IJ (t^m)" e {t*m) 

xeM xeM 
Denoting generalized velocities by 



n+l 



,• dx'' ,• de' 

V = -rr, v'a = 



a 



dt dt 
and their conjugate canonical momenta by pi, p^i, we have that 

py+p^iV^A = PiV + P^iV^A, 

and 

V' = v\ V'a = v'a + rjk{x)e^AV^, 

Pi=Pi- e^AP^k^'^oU P^i = P^i- 

Let us observe the following anti-dualism: it is easily seen that and P"^, are 
connection-independent, whereas Pj and V^a depend explicitly on F. Geometric 
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reasons for this are that that V is simply the 7r-projection of gcnerahzcd velocity 
Q{t) to M and P^j are components of the linear functional on Tg(^^)Fx(^^)M 
{x{t) = 7r(g(t))) which is simply the usual restriction of some functional on 
Tgi^qFM to the subspacc Tg^^-^F^^^-^M . The restriction procedure is obviously 
connection- independent . 

Remark: v^a are not components of vectors in T^^M, and similarly, pi are not 
elements of covcctors. Only and p'^i arc so. The total systems [v'^,v^a), 
{PiiP^i) represent the higher-floor vectors and covectors in T^FM. The con- 
nection-dependent systems a, Pi are respectively components of vectors and 
covcctors in T^M (elements of F^M itself and F*M). 

Quite a similar reasoning works for infinitesimal gyroscope, i.e., for the bun- 
dle F{M,g). The situation there is more complicated technically, because the 
quantities e^A are no longer independent, i.e., they satisfy the orthonormality 
conditions 

Qije^e^B = Sab- 

This fact creates some problems in equations of motion; to overcome them one 

uses an auxiliary technical tool, namely, a field of linear orthonormal frames 
defined all over the manifold M. We return to this question later on. 

In mechanics of extended afRnely- and metrically-rigid bodies one uses the 
concepts of affinc velocity (Eringcn's "gyration") and affinc spin (affine momen- 
tum, hypermomentum) . They may be as well defined for infinitesimal objects, 
i.e., for internal degrees of freedom. 

Every tensor object in the tangent space T^FM or TeF{M,g) may be ex- 
pressed in terms of its components with respect to the frame e itself. In partic- 
ular, this concerns the translational velocity V and the system of internal veloc- 
ities (. . . , Va, ■ ■ •)• Similarly, the covariant conjugate momenta P, (. . . , P^, . . .) 
may be expanded with respect to the dual co-frame e. So, we have that 

V = V^eA, Va = esQ^A, 

or explicitly, using differentiation symbols, 

^ \ ' Dt / Dt 

This expression is exactly the co-moving afHne velocity known from the theory 
of extended afBne bodies in a flat space. V"^ are co-moving components of the 
translational velocity. Geometrically, these objects belong to the matrix spaces, 
i.e., 

fl e L(n, M) ~ GL(n, R)', F e M" ~ L(l, n; R). 

But having e instantaneously fixed, we can associate with these numerical ob- 
jects the corresponding quantities in the tangent space T^M, where x = 7r(e) 
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(e S FxM). For V this is just the usual translational velocity V € T^-M: 

V' = — 
dt ■ 

The quantity 51 and the frame e give rise to the object 

i.e., explicitly in terms of derivatives: 

o DeA De^A a z f]A b 

n = -^(S,e , n^^—^e j^CA^ Be j. 

Geometrically, they are linear transformations in tangent spaces, i.e., 
n e L{T^M) ~ T^M (E) T*M = GL(r^Af)', 

where, obviously, x — 7r(e) (e € FxM). 

Just as in mechanics of extended affine bodies, fl^j are spatial (laboratory) 
components of the affine velocity, i.e., they represent what Eringen used to call 
"gyration" [101 [11]. Geometrically it is important that fl, belong respectively 
to Lie algebras of the groups GL(n, M), GL{TxM). 

The same may be done for the covariant canonical momenta P and P"^: 

P^Pac^, P^^t^Be^, 

i.e., 

Pa - {P, ca) = P^e'A, S^b = (P^, es> - P^^e'B. 

Just as previously, P G M", E E L(n,M), but strictly speaking one means here 
R" as identified in a Cartesian way with its own dual M"*; similarly, L(n,R) 
plays here a role of the Lie co-algebra L(n, R)* of GL(n, R). One does not notice 
this subtle distinction because L(n,M)* and L(n,R) are canonically identified 
via the pairing 

{A,B) = Tt{AB). 

In analogy to mechanics of extended bodies in flat spaces we say that the 
matrix elements of S are co-moving components of the afhne spin (hypermo- 
mentum). Just as in flat-space theory, the laboratory (spatial) description is 
based on the quantity 

Es^A ^ B i T-tA i \^A B 

= 2j BCA'^e , 1j j = e aP j = e a^ Be j. 

So, S'^B are Hamiltonian generators of the structural group GL(n,M) (micro- 
material global transformations) ([T|), ([7]). And similarly, are Hamiltonian 
generators of the local group of microspatial transformations (jl]), ([5]), (O. 

If the configuration e G FM is fixed, then the velocities , V'a contain 
exactly the same information as F*, f2*j or V^, ^I^b- Similarly, the momenta 
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Pi, P^i are equivalent to Pi, E*j and Pa, 'S^b- Using more sophisticated terms 
we would say that the state space TFM is naturally diffeomorphic with the 
manifolds 

TnFM := \J © T.M ® L(r,M), 

Similarly, the phase space T* FM is canonically diffeomorphic with 
T*FM y F^M ®T*M ®HT^M), 

T£FM := y F^M®R"®L(n,M). 

xeM 

Let us notice that there is a subtle distinction between T^FM and T^FM 
indistinguishable in our simplified notation. Namely, in the latter space the 
R"-term consists of row numerical vectors, i.e., linear functionals on the space 
of numerical column vectors (so, strictly speaking, one deals with R"*). And 
L(n,R) is canonically isomorphic with L(n,R)*. (Remark: this isomorphism 
is canonical for L{V), L{V)* , where V is an arbitrary linear space. This is no 
longer the case for V, V* themselves.) Obviously, the mentioned identifications 
between state manifolds may be interpreted as natural ones only on the basis 
of some fixed affine connection F on M. 

The cotangent bundle P = T*FM carries a natural symplectic structure 
[U 1331 [55] and is used as the mechanical phase space of our problem. As just 
mentioned, the affine connection T fixes some diffeomorphism of P onto the 
manifold V ~ T*FM . This diffeomorphism enables one to carry over to V the 
intrinsic symplectic geometry of P = T*FM . But there is one delicate point 
which, when overlooked, may lead to serious mistakes. Namely, unlike the pre- 
viously used coordinates in T*FM, the quantities x% e^A, Pi, P^i fail to be 
canonical (Darboux) coordinates for the symplectic structure F-transferred to 
T*FM . After some essentially easy although sometimes technically embarrass- 
ing calculations one obtains the system of basic Poisson brackets. So, just like 
in mechanics of extended affine bodies, we have the obvious rules 

{x\x^ = (), {e'A,e^B} = 0, {x\e^A}^Q, {P^,P^,} = 0, 
{P^,a;^} = 0, {e^A,P^,} = ,5^,<5^^, {x\Pj} = 5',. 

However, further on one obtains more complicated expressions explicitly depen- 
dent on the F-geometry of A/: 

{P„ P,} = S^7^'HJ, {P^, P^} = -P-^fcr^., {P„e^A} = e^F^H, 
where, obviously, IZ denotes the curvature tensor of F; we use the convention 

hij — J- bj.i hi.j ~r 1 ci^ bj cj^ hi 
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(comma, as usual in the tensor calculus, denotes the partial differentiation with 
respect to the indicated coordinate). The latter three Poisson brackets mean, 
roughly speaking, that Pi are Hamiltonian generators of parallel transports of 
our state variables. One can also show that 

{I]%-,I]M = (5'zS]^--,5'=,I]^. 

In these formulas we recognize structure constants of the linear group. As 
expected, this means that E^^ are basic Hamiltonian generators of (|4|), i.e., 
roughly speaking, of the system of groups GL(Tj;M). It is easy to obtain the 
relationship 

that is also compatible with the mentioned interpretation of Pi as Hamiltonian 
generators of parallel transports. 

If some function F depends only on the configuration, i.e., on the FM- 
variables (but not on Pi and P^i)^ then 

{^^,,F}^^E^,F^-e^.§-, 

where 

F j = e KS jF L = e kt——- 
The Poisson brackets involving co- moving components are as follows: 

{Pa, Pb} = ^'^lH'^kab - 2PkS'^ ab, {^^b,Pc} - -PbS^c, 

{Y.^ B,^'" d} = B^^D — ^^d'Y'" B, {^^Jt'^^b} = 0, 

where TZ^ kab, ab are respectively co-moving components of the curvature 
and torsion tensors of F with respect to the instantaneous internal configuration 
e, i.e., 

/<• KAB ■— e iK jmnd KG A^ B, 

cK K ni j m ci ^ f T^i F'* ^ 

AB ■— c iO jjyiti A^ Bi 'J jm — 2 ■'™ " 

Let us observe a characteristic difference between the spatial and co-moving 
representations of Poisson brackets. Namely, the latter ones are "almost" iden- 
tical with the basic commutation relations (structure constants) for the affine 
group GAf(n,K). The "almost" concerns the brackets {Pa,Pb} which do not 
vanish if the connection F is not completely flat (in the sense that both the 
curvature and torsion tensors vanish). The co- moving brackets enable one to 
interpret Pa as a Hamiltonian generator of parallel transports along the A-th 
legs of the frames e. 

Let us also note certain additional and convenient Poisson brackets. If F 
depends only on the configuration, then 

{S^B,F} = -i?%P = -e*B^, 

de*A 

{Pa.F) = -HaF, {P„F} = -H,F, 
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Again we conclude that Pi arc Hamiltonian generators of parallel transports 
along the i-th coordinate axes, and Pa are generators of parallel transports 
along the ^-th legs of the frames e. 

In the theory of an extended affinely-rigid body in Euclidean or affine space 
the quantities Pi, P^i, S'^b have a very natural geometric interpretation 
based on some transformation groups acting in the configuration space. Let us 
remind that in the flat affine case the configuration space FM trivializes to the 
Cartesian product M x F{y), where F{V) is the manifold of frames in V. The 
most important transformation groups acting in M x FiV) are the following 
ones: 

(i) spatial translations (x, e) {ty{x)^e); they shift the centre of mass 
X £ M along the vector v gV without affecting the internal configuration 
e€ F{V). Analytically: 

{x\e'A) ^ {x' +v\e'A) ■ 

(ii) additive translations of internal degrees of freedom: any ^ S V" 

gives rise to the mapping {x, e) i-> (a;, e + ^), i.e., 

{x\e\)^{x\e\ + eA). 

Of course, such transformations act only locally in M x F{V), because 
they may produce linearly dependent n-tuples of vectors from independent 
ones. The Hamiltonian generators are given by P^i. 

{Hi) spatial £iffine transformations of internal degrees of freedom: any 
L e GL(y) acts on the configuration space as follows: 

(a;; . . . ,eA, ■ • ■) 1-^ (x; . . .,LeA, ■ . .) 
or shortly (x, e) {x,Le). Analytically, 

{x\e'A) ^ {x\Uje^A)- 

The centre-of-mass position is not effected. The Hamiltonian generators 
of this group are given by Because of this the object is referred 

to as affine spin. Ii g G V* ^ V* is the metric tensor of V, then the 
quantity 

S^j = S*j — gjk^^ig^^ ^ 

i.e., the doubled (/-skew-symmetric part of E, is the usual canonical spin 
generating rigid spatial rotations of internal (relative) degrees of freedom. 
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material afRne transformations of internal degrees of freedom: 
any not singular matrix L G GL(n, R) acts on the configuration space as 
follows: 

{x;...,eA,---)'-^ (x;..., cbL^a, • • •)• 

As in (iii), we use the shorthand (x, e) {x,eL). Important in this 
context: do not confuse conceptually linear mappings in V with matrices! 
This group is generated in the Hamiltonian sense by S'^b, i.e., by co- 
moving components of the hypermomentum. Let us observe that the above 
right-acting transformations may be interpreted as ones induced by Unear 
mappings L e GL(n,M) acting in the "material space" R". As usual, the 
frames themselves may be interpreted as linear mappings e : M" V. The 
label space M" is endowed with the standard metric S, sometimes denoted 
also by rj, to stress the link with the usual formulation of the mechanics 
of extended afRne systems. The (5-skew-symmetric part of T,^b, 

V^B '■= — Sbc^*^ d5^^, 

i.e., vorticity, generates rigid material rotations of the body. 

total afRne transformations in space. They act both on the transla- 
tional and internal degrees of freedom. Let (p : M M he an arbitrary 
afEne transformation acting in M, and L[ip] : V ^ V he its linear (ho- 
mogeneous) part. Explicitly, for any pair of points a,b G M, the radius 
vector (fi{a)(fi{b) is obtained from the vector ab through the action of L[ip]: 

ip{a)ip{b) = L[<f]ab. 

The action of (p on the configuration space M x F{V) is given by (a;, e) i— > 
{(p{x), L[(p]e). Analytically, in Cartesian coordinates x', 

(x'.c'a) ^ {L'jX^ +a\Vje^A) ■ 

Infinitesimal Hamiltonian generators arc given by i^Pi^J'^j)^ where the 
total canonical affine momentum J^"^ j with respect to the origin of coordi- 
nates is defined as follows: 

i7 j = a; + S j , 

i.e., it consists of the "orbital" and "internal" parts, respectively, x' Pj 
and E*j. Taking the skew-symmetric part of the above expression we 
obtain the usual splitting of the angular momentum onto its "orbital" and 
"internal" (spin) parts: 

= Uj + S^j, 

where, obviously, 

^ j ~ ^ j 9 9ji'^ k — \J j 5 

= x'P,-g'''g,ix'Pk=x'Pj^xjP\ 
D j — Zj j — g yjl^ k — ^ j ~ ■ 
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(m) total afRne transformations in the material space. If the configu- 
ration space of extended affinely-rigid body is identified with M x F{V), 
then, as mentioned, the material space may be identified simply with M". 

The current position of the constituent labelled by a G R" is displaced with 
respect to the instantaneous position of the centre of mass by the spatial 
vector u with coordinates = eVct^. The afEne group GAf(n, R) in M" 
is canonically identical with the semidirect product GL(n,R) R", i.e., 
with the set of pairs {B, b) multiplied (composed) according to the group 
rule 

{BiM){B2M) = {BiB2,bi + Bib2). 
The action of {B, b) on the configuration {x, e) is given by 

(a;,e) ^ {teb{x),eB), 

where in the expression eb, e G F{y) is identified with a linear isomor- 
phism of R" onto V. Analytically, 



The corresponding system of infinitesimal Hamiltonian generators is given 



If (M, r) is not flat, then the above picture of the transformation groups 
changes in an essential way. There is no counterpart of the n{n + l)-dimensional 
group of affine transformations acting in Ad . In particular, there is no concept 
of spatial translations and radius vectors. Only the infinite-dimensional group 
Diff(M) of all diffeomorphisms of M onto M is well defined. As a rule, its 
elements do not preserve the parallel transport and covariant differentiation. 
There are no extended affinely-rigid bodies and affine degrees of freedom may 
be considered only as internal ones. Translations in the micromaterial space 
R" are well defined, however, they lose the physicial meaning that they had in 
the theory of extended affine bodies. The only transformations which survive 
are those dealing with internal degrees of freedom only, without any affecting 
of translational motion in M. Let us describe them in some details: 

(i) additive translations of internal degrees of freedom. They act in 
any fibre F^M = tt~^{x) C FM of the bundle of frames exactly as their 
flat-space counterparts do in F{V). There is, however, an essential novelty 
when these transformations are considered globally in FM , not in separate 
fibres. Namely, in a curved manifold there is no canonical identification 
between different tangent spaces. Therefore, the concept of "the same" 
translation in different fibres is missing, and the group becomes infinite- 
dimensional. Every ordered n-tuple of vector fields in M, 




by 



M B X ^ ax) = {■ ■ ■ ■ ■ ■) & {T,M) 



n 



gives rise to the local transformation of FM 



F^M B e>-^ e + ^(x) e F^M 



(8) 
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(local, because some frames are mapped into linearly dependent n-tuples 
of vectors). Obviously, the sum in the last formula is meant pointwisely, 
i.e., for any A — l,n the element ba is replaced by ca + £,a{x)- Canonical 
momenta P^i are Hamiltonian generators of such transformations. More 
precisely, for any ordered n-tuple of covector fields Xa on M, the quantity 
X^aP^i generates some one-parameter group of transformations ([5]). 

(m) spatial affine transformations of internal degrees of freedom. Here 
we are dealing with the same phenomenon as previously, namely, the n^- 
dimensional Lie groups GL(ra;M) acting separately in tangent spaces T^M 
give rise to some infinite-dimensional group of transformations acting glob- 
ally in FM . This group is "parameterized" by mixed second-order tensor 
fields T on M, cf. (HI, ([5|), ([6|). The quantities 5]*^ are Hamiltonian gen- 
erators of this group. More precisely, for any mixed second-order tensor 
field a on M, the quantity ^Y? i generates some one-parameter group of 
transformations (|4]), ([5]), ([6]). There is no escape from the infinite dimen- 
sion because in a curved space (M , F) there is no canonical identification 
of tangent spaces at different points. 

iiii) micromaterial afRne transformations of internal degrees of free- 
dom. Now, in a complete analogy to the extended model in a flat space, 
the n^-dimensional Lie group GL(n, K) acts on the configuration space 
FM according to the usual rule ([1]). The quantities S'^b are Hamiltonian 
generators of the corresponding Lie groups of extended point transforma- 
tions. 

One should stress that the quantities b describing micromaterial trans- 
formations may be defined as constants, but they need not be so. Namely, 
following the pattern developed in gauge theories we can replace the matri- 
ces L by GL(n, M)-valued functions on M. They act on the configuration 
space FM according to the rule ([3]), ([7]). Any matrix- valued function 
a : M — > L(n, R) gives rise to the one-parameter group of transforma- 
tions with the Hamiltonian generator ol^b^^ A- In this way we obtain 
again, just as in previous examples, some infinite-dimensional functionally- 
" parameterized" transformation group. 

We shall almost not deal with the diffeomorphism group Diff(Af) or its spe- 
cial (volume-preserving) subgroup SDiff(Af) as a symmetry of our model. It 
does not play any essential role in mechanics of material points with internal 
degrees of freedom. On the contrary, it is rather relevant for the theory of 
micromorphic continua. 

Here let us only mention that any diffeomorphism / : M — > M acts on the 
configuration space FM according to the rule 



F^M 3 Df^ ■ e 



i.e., analytically, 
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where the functions /'(x) provide analytical description of / (they express co- 
ordinates of /(x) in terms of those of x). 

Due to the essentially local ( "x-dependent" ) character of transformations 
(HI), (HI, ([7]), dHl), and admissibly local character of their action on velocities 
and canonical momenta is (respectively, may be) different from that could be 
naively expected on the basis of analogy with the flat-space theory. 

Additive translations in internal degrees of freedom ^ affect internal veloc- 
ities as follows: 

where, as previously, = dx^ /dt denotes the translational velocity in M. This 
rule reduces to the invariance of (. . . , V^i, . . .) (characteristic for the flat-space 
theory) only when all vector fields are parallel in the F-sense. Translational 
velocity is invariant under transformations of internal degrees of freedom, 'V" = 
. On the contrary, the covariant canonical momenta transform according to 
the rules 

Microspatial affine transformations of internal degrees of freedom (jl]), (O, (H]) 
also do not affect translational velocity, i.e., 

'V = V\ 

but the internal velocities are transformed in the following way: 

'V'a = T'jV^A + (Vfer,) e^A. (9) 

This implies that 

'n'j = r,r!'„T"i'"j + y'^ (VfcTV) p-'^j 

= r,r!'„T'i™j + (VyP%„)T-i™j, (10) 

and similarly, 

'^^B - ^^B + e^iT-'', i^vT'j) e^B- (11) 
On the contrary, translational canonical momenta Pi suffer the transformation 

'P, = P, - S^T^i', V.T^fc, (12) 

whereas the internal ones obey the well-known global rule: 

'pA _ pA rp-lj 

therefore, 

f\^i rpi x^k rr!—lm /\^A \^A 

^ j ~ J- m-L j, ^ B ~ ^ B- 

The formulas pop . (fTTj) . (fT^ reduce to the corresponding global rules from 
the mechanics of extended affine bodies if and only if the field T is F-parallel, 
i.e., 

Vfe^j = 0. 
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The global micromaterial transformations ([T]) act just like in the flat-space the- 
ory of extended affine bodies: 



= V\ 'V'a = VbL^'a. 



11 

therefore, 

= ri*j, '^^B — L ^^c^'^ dL^ Bi 

^ j — ^ ji ^ B — B- 

For local, i.e., x-dependent micromaterial transformations we have that 

'^'j = n',+e:B{L^A,kL-^^c)e^jV\ 

where comma denotes the partial differentiation. Again these formulas reduce 
to the preceding ones when the field L is constant, i.e., micromaterial transfor- 
mations are global. By duality rule, 

'P,'V' + 't'^A'^^B = P^V' + t'^A^^Bi 

one can easily show that 



'Pi — Pi — ^^KL^B.iP , 



Ki^ BA^ 
B ~ D^J Bi 

'P^, = L-'^bP"" 



3 Metrical concepts 

We have described above kinematics and canonical formalism for afRne model 
of internal degrees of freedom in a manifold M endowed with a general affine 
connection T. The metrical concepts were yet not used at all. They become 
necessary when we aim at constructing dynamical models. And of course, in 
realistic physical problems one usually deals with some metric structure. For ex- 
ample, in the standard General Relativity, i.e., in the Einstein theory of gravita- 
tion, the dynamical metric tensor of the four-dimensional space-time is used for 
describing gravitational field. Even in more sophisticated alternative theories, 
like gauge models of gravitation, there exists always some physically relevant 
metrical aspect of the gravitational field. 

So, from now on, we shall usually (nevertheless, not always) assume that the 
manifold M is endowed both with the metrical tensor g and affine connection 
r. In principle, one can consider the most general structure, where F and g are 
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unrelated, independent on each other. Nevertheless, both from the geometrical 
and physical point of view of particular interest are situations when some kind 
of compatibility between afSne and metrical structures is assumed. 

Riemann-Cartan space is one in which all metrical relations like length of 
vectors, angles (in particular, orthogonality) are preserved by parallel trans- 
ports, i.e., 

Vfefifij = 0; 

the covariant derivative of the metric tensor vanishes. It is well known from 
differential geometry that then 

r%fe = \ I + K''jk = I \ + S'jk + Sjk + Skj\ (13) 



where 



1 



jj^ ( — {9mj,k 9mk,j gjk,m) 



is the Levi-Civita (Christoffel) symbol, 

S'jk = T\jk] = 2 i^'jk - r'fcj) 

is the torsion tensor, and all indices are moved from their natural positions 
with the help of the metric tensor g. The quantity K^jk is referred to as the 
contortion tensor; it satisfies 

K'jk + K/k = 0. 

The structure {M,r,g) reduces to the Riemann space when it is torsion- 
free, S^jk — 0, i.e., the object T is symmetric and automatically coincides with 
Levi-Civita symbol, 

i 



jk 



jk 



In this model parallelism is derived from the metric concepts. 

Riemann-Cartan- Weyl space is one in which angles between vectors, but not 
necessarily their lengths are preserved, i.e., 

Vfeffij = —Qkgij, 

where Qk is referred to as the Weyl covector, and its contravariant counterpart 
Qk _ gkjQ^ as the Weyl vector. It is known from differential geometry that in 
Riemann-Cartan- Weyl spaces 



where ^ 

K^^jk = S''jk + Sjk + Skj^ + - {5^jQk + ^^Qj — gjkQ^) 
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with the same as previously convention concerning the raising and lowering of 
indices. 

If r is symmetric, we are dealing with the Weyl space, 

r^-fe = {jk] + l (^'^-^fc + ^'"^^ - 9jkQ') ■ 

In the micromaterial space M" the standard Kronecker metric Sab is used. If 
for any reason wc admit non-orthogonal rcctiUncar coordinates in M", we shall 
use the symbol tjab to denote the micromaterial metric. 

For any e € F^M the corresponding Green deformation tensor in R" is given 

by 

G[e]AB = gij{x)e\e^ B 
or, using the obvious abbreviations, 

G[e] =e* -ge W* ®W* . 

Similarly, the Cauchy deformation tensor C[e] &T*M ®T*M is given by 

C[e]ij = "nABe^iB^ j = dABe'^iC^j, 

or symbolically, 

C[e] = g* • ?7. 

The corresponding contravariant inverses will be denoted by G[e], C[e], 

AC P^ik /~i 

LrCB = B, (-^kj = j. 

Obviously, 

G[e]^^ = e^ie'^jQixf^ C[e] = e^^e^B??^^ = e'^e^s^^^. 
Lagrange and Euler deformation tensors are respectively given by 

S[e] = \{G-ii), e[e\ = \{g-C). 

Obviously, the quantities Gfej^s are scalar products of vectors ba, sb- 
G[e]AB = g{x) {eA,eB) = (e^les) = G[e]BA- 

We say that motion is metrically rigid if £ (equivalently e) vanishes along all 
admissible trajectories. In other words, only such configurations are admissible 
that 

G[e]AB = Vab = Sab- 
Equivalently, this means that 

C[e] = g{x), e € F,M. 
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Infinitesimal affinely- rigid body becomes then the infinitesimal gyroscope, i.e., 
the configuration space FM becomes restricted to the constraints submanifold 
F{M, g) C FM consisting of (;-orthonormal frames. More precisely, one ad- 
mits the connected submanifold F^ (M, g) consisting of ij-orthonormal frames 
positively oriented with respect to some fixed orientation in M . 

As mentioned, constraints of gyroscopic motion may be described analyti- 
cally by any of the following systems of equations: 

gijc'-Ae^B = VAB = Sab, (14) 

VABe^iS^j = ^ABe^^ie^] = 9ij, (15) 

obviously, these systems are equivalent. 

Subjecting the first system to the operation D/Dt, i.e., covariant differenti- 
ation along the orbital trajectory, and contracting the resulting equation with 
e^^iC^j, we obtain that 

9M^"i + gu^'k = -V'^Vmgki = -^vgu- (16) 

Equivalently, differentiating the second system in the D / Dt-sense, we obtain 
the equivalent form 

VAC^'^B + VBC^'^A = - (V^^V^g.j) e^e^B = - (Vy^.j) e\e^B- (17) 

Let us notice that the ZJ/ZJt-differentiation of the first system is identical with 
the usual d/dt-differentiation, because from the point of view of geometry of 
AI the left- and right-hand sides are scalar quantities (they are tensors in the 
micromaterial sense of R"). 

If r is metrical, i.e., if (M,T,g) is a Riemann-Cartan space, in particular, 
just a Riemann space {M,{},g), then the right-hand sides of equations ([TB|) . 
(jl7p vanish and f2, Q become respectively g- and 77-skew-symmetric, 

gkm^"'! + gim^'^k = 0, vkm^"l + vlm^^'k = 0. 

If the standard coordinates in M" are used, tiab = 5ab, then the last condition 
means simply that f2 is literally skew-symmetric. 

Therefore, il, Q. are respectively elements of Lie algebras SO{TxM^gx)' , 
SO(n, R)' of the corresponding orthogonal groups. They are gyroscopic angular 
velocities respectively in the spatial and co-moving (material) representations. 
For a unconstrained affine motion, when f2, are general linear mappings, an- 
gular velocity may be defined as the corresponding g- or 77-skew-symmetric part 
of the affine velocity, i.e., 

^'3 ■■= li^'.-g-'^gji^'k), 

In the case of rigid motion they are identical respectively with r2*_,- and VL'^b- 
But attention! The general rule 
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is not valid any longer for tj'j, lo^b, i-e., 

unless the motion is metrically rigid (gyroscopic), i.e., (|14p . (fT5|) hold. Therefore, 
in non-gyroscopic motion lo^b fail to be co-moving components of the angular 
velocity. The same concerns kinematical distortions, 

d^B l{n^B+V^^VBDn"c). 

As seen from the formulas ([TS]). p?)) . if affine connection F used in the 
definitions of fl and is not metrical (i.e., (M, F, g) is not the Riemann-Cartan 
space), then and fl are not skew-symmetric in the g- and 77-sense (i.e., are not 
elements of the Lie algebras of SO{TxM, g^), SO{n, M)) even if motion is purely 
gyroscopic. Therefore, it is confusing to interpret them as angular velocities. 
Taking their g- and ?7-skew-symmetric parts also does not seem convincing. If 
F is metrical, i.e., — 0, then, as mentioned, f2 and D, are respectively g- 
and r/-skew-symmetric. Nevertheless, their particular form depends explicitly 
on the torsion tensor S. So, no doubt, the use of gyroscopic concepts is clean 
only when the Levi-Civita afhne connection {} is used. 

One can wonder whether gyroscopic motion could not be defined without 
any use of the fixed metric tensor g on M, thus, without any problems like 
those mentioned above (definition of angular velocity, T-g compatibility, and 
so on). Apparently, the idea might look both not physical and mathematically 
inconsistent. However, although the physical usefulness question is still open, 
the mathematical correctness may be easily shown. In a sense, gyroscopic con- 
straints may be defined on the basis of afhne connection structure (M, F), quite 
independently of the metric tensor concept. However, the micromaterial metric 
VAB (Sab) may be used. Namely, it gives rise to the Cauchy deformation tensor 

Cij — fjAB^^i^^ j 

and its inverse ^ 

^e^AC^Bv"^^- 

Obviously, C[e] G T*^^^M®T*^^^-^M has all formal properties of the metric tensor 
in the instantaneous tangent space T7r(e)Af (symmetry and positive definiteness) . 
But it is defined only at the point x = 7r(e) € M and is not induced by any 
metric field g living globally all over in M . And without such a field even the 
very term "deformation tensor" is not very adequate, because we do not have 
any metrical standard which might be compared with C[e]; thus, we cannot 
decide to which extent e G FM "deforms" 77. Nevertheless, it is meaningful to 
say that some motion in FM is free of deformations when C[e] is covariantly 
constant along the curve describing translational motion in M, i.e., 
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After simple calculations we obtain the mutually equivalent conditions 




(18) 
(19) 



i.e., VL and are skew-symmetric respectively with respect to ry and C used 
as metrics. There is however an important point, namely, if (M, F) is non- 
Euclidean (in the local sense), then the above constraints are aholonomic. If 
(M, r) is locally Euclidean, these constraints become quasi- holonomic, i.e., the 
manifold FM is foliated by a family of n(n-f l)/2-dimensional integral manifolds 
of the above Pfaff systems. Any of these mutually disjoint manifolds is a possible 
configuration space of rigid motions. And any particular choice is equivalent to 
fixing some metric tensor field g on M . And indeed, in an n-dimensional linear 
space the manifold of possible metric tensors is n{n -\- l)/2-dimensional, thus, 
leaving n{n — l)/2 independent parameters in the manifold of linear frames 
reducing the metric to its standard Kronecker-delta form. Together with n 
translational degrees of freedom we obtain exactly n(n-f l)/2 degrees of freedom 
of a rigid body moving in n-dimensional Euclidean space. 

It is quite a different question whether the Cauchy tensor C[e] may be mean- 
ingfully used as a metric-like tensor of the instantaneous tangent space T^(e)M. 
We return later on to the question of possible physical applications. In a mo- 
ment we stress interesting geometrical aspects of aholonomic constraints ^8]) . 



4 Dynamical afRne models 

Let us now turn to describing dynamical models. We begin from the models 
of kinetic energy, i.e., roughly speaking, metrics (Riemannian structures) on 
the configuration space FM . When we deal with internal degrees of freedom, 
the problem becomes delicate because the d'Alembert method of deriving the 
kinetic energy from the model of extended affinely constrained system becomes 
unjustified, not reliable, perhaps just misleading. It is so even in models of 
essentially internal degrees of freedom in a flat space, but in curved manifolds 
the problem becomes very essential, fundamentally embarrassing, and the only 
reliable method is one based on appropriate symmetry principle. As we saw in 
|49| , it is the case also in some non-standard applications in the usual continuum 
mechanics and the dynamics of structured bodies and defects. 

Nevertheless, to begin with, we discuss as first the models (metrics on FM) 
following formally the expressions known from the d'Alembert mechanics of 
extended affine bodies in flat spaces, and in any case based on simple analogies. 
Later on we discuss mathematical structure of more general models following 
the non-d'Alembert ideas in a flat space [15]. And finally, some perspectives of 
physical applications will be suggested and preliminarily discussed. 

Repeating formally the d'Alembert expression from mechanics of affine bod- 
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ies in flat spaces, we obtain that 



2^^ dt dt 2^^ Dt Dt ^ ' 

In this formula the descriptors "tr" and "int" refer obviously to the translational 
and internal parts, m > is the mass, and the symmetric and positively definite 
micromaterial tensor J e M" (g) R" describes the internal inertia. The main 
difference in comparison with the flat-space formula for extended bodies is that 
now the covariant along-curve derivatives of "directors" ca are used. And there 
is of course one subtle point; namely, now J is not derived from the model 
of constrained extended system, but just postulated as something primary. In 
certain calculations the following equivalent expressions are convenient; 

Ttr = yG[e]^By^i>^, (21) 

Tint = ^G[e]KL^''An'^BJ^''. (22) 

Obviously, now the coefficients in quadratic forms of co-moving velocities depend 
on the internal configuration e. Alternatively, Tint may be expressed in spatial 
terms: 

Tint = ^g,,n^n^iJ[ef', (23) 

where J[e\ is the configuration-dependent spatial representation of the internal 
inertia: 

jt^]kl _ J jAB 

J[e\ — e AC bJ 

This second-order moment (quadrupole of the internal inertia) in extended body 
dynamics is often used in nuclear physics. 

Remark: Tint in (EOl) , (EH) is invariant under translations ^ if and only if the 
field ^ of additive translations is parallel under F, i.e., Vp]^ = 0. 

Denoting the system of generalized coordinates of affine bodies by 

(...,g^,...) = {...,x\...;...,e'A,---), 
and writing symbolically (|20p in the following form: 



T = -Q^I.v 



1 dq^' dq" 



we see that the underlying Riemannian metric Q on F{M) is fiat if (M, g) is 
locally Euclidean. It is no longer the case in a curved manifold. 

In non-dissipative models with the velocity-independent Lagrangians L = 
T — V{x,e) the resulting Euler-Lagrange equations 

dtdqi^ dqf ^ ' 
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may be written down in the following form: 
DV 



DV"- 



-Krr^^al^^E^jAB ^pa (25) 

Dt Dt ^ ' 

eV^^^J^^ = -6^5*""^^^'^^ + A^"', (26) 

Dt^ ^ Dt Dt ^ ' 

where the meaning of symbols is as follows: 

ic^a -pa 



a 
be 



abfdV , dV 

dV 



and obviously, all shifts of tensor indices from their natural positions are meant 
in the sense of g. 

Let us stress an important fact that in general the translational force F" 
docs not equal —g'^^'dV/dx''. The equality holds only when V does not depend 
on internal degrees of freedom, i.e., when y is a 7r-pull-back of some function 
on M, or if connection T is flat and local Cartesian coordinates on M are used. 
If V is not a 7r-vcrtical function on M, then dV/dx^ is not a covariant vector in 
M at all, and g'^^dV/dx'^ fails to be a contravariant M- vector. Let us observe 
that F may be written as follows: 



dV 
da* 

The CO- moving representations oi F, N are given as follows: 
Fa = Fie\ = -HAV, 

oe^A 

Remark: the contravariant objects F^ = eAiF"^, N^^ — e'^ie^jN'^^ have the 
forms 

i.e., the lower-case indices are raised with the help of the Green deformation 
tensor. 

When non-potential interactions, e.g., dissipative ones, are admitted, the 
above general form of equations of motion in principle remains valid, however 
the dynamical terms f , TV"'' must be replaced by more general expressions pos- 
tulated on independent basis (e.g., additional friction forces linear in generalized 
velocities with symmetric negatively-definite coefficients matrices). 
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It would be technically very difficult to derive equations (pS)) directly from 
the second-kind Lagrange equations ([24]) . It is much more convenient to use the 
canonical formalism and Poisson-bracket techniques. For the potential systems 
with Lagrangians L ^ T ~ V (a;*, c^a), the Legendre transformation 

dL dT 

leads to the Hamiltonian 



where the kinetic term 



H = T + V (x\e 



has the form 



T = ^9''m + IjABP^^P'^jg''. (27) 



Let us remind that J is the inverse of J, 

J AC 7 

J JCB — O B, 

and the explicit expression for Legendre transformation reads that 
FIT r)T 

The formerly quoted basic Poisson brackets enable one to write down explicitly 
Hamiltonian equations of motion, 

rIP rlP^- 

which after some manipulations may be reduced to the form (j25p . Let us stress 
that the covariant derivative in (j25p is meant in the sense of F, not in the g-hevi- 
Civita sense; an important fact to be kept in mind when they do not coincide 
(i.e., when {M,T,g) is not a Riemann space). 

As expected, equations of motion (I25p simplify in a remarkable way when 
(M, F, g) is a Riemann- Cart an space, = 0. In this case K, becomes so-called 
contortion, 

be ^ ^ be + ^bc + Ocb , 

and equations (^5)) reduce to 

nya 1 

m-^ - -S^n^k^jV^ +2mV''V'Sbc'' + F'', (28) 
e-^K^^^J^^ = N^'. (29) 
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In the first (translational) equation only the g-skew-symmetric part of i.e., 
S''j , survives, because in Riemann-Cartan spaces the curvature tensor is g-skew- 
symmetric in the first pair of indices. Besides of the usual external force F°', 
the right-hand side of involves two geometric forces describing the coupling 
between spatial geometry and kinematical quantities of the particle motion. 
Namely, translational velocity is quadratically coupled to the torsion, and spin 
is coupled to the curvature. This is in a nice way compatible with the geometric 
interpretation of curvature and torsion respectively in terms of rotations and 
translations. The geometric force 

is so-to-speak magnetic-like in the sense that due to the corresponding anti- 
symmetries of R and S they are g-orthogonal to velocities and do not do any 
work, i.e., 

n n 

yab^ ^ gcom " i 

therefore, they do not contribute to the energy balance. In particular, they do 
not accelerate the absolute value of F", || — \/gabV°-V'', but only change the 
direction of V"", i.e., give rise to the bending of translational trajectory. Even 
in the force-free case the motion is not geodetic; there is some link between this 
phenomenon and geodetic deviation. Let us observe that (|28|) may be further 
simplified to the form 

Dt 2 

where -D[g] denotes the along-curve covariant differentiation in the Levi-Civita 
sense. However, if (M, F, g) is non-Riemannian, i.e., the torsion does not vanish, 
this simplification is a rather seeming one, because the curvature tensor TZ of 
r contains both the term corresponding to the Riemannian curvature of {} 
and the terms involving torsion. And in the internal equation ((29|) there is no 
simplification at all. Obviously, in the purely Riemannian case, when the torsion 
does vanish, we obtain the maximally simple and clear system of equations of 
motion: 

i5^7^'fc^l^^ +F^ (30) 
7v°^ (31) 

and now both the covariant derivative and the curvature tensor are meant in 
the (7-Levi-Civita sense. 

It is instructive and convenient to write down these equations in terms of 
some balance laws. Let 

dt ' 

T^ab a T^b jKL a L jKL 



771- 



Dt 
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denote respectively the kinematical linear momentum and kincmatical affine 
spin (kinematical hypermomentum) , just as in [4^. Then ((30l) . (|3T|) may be 
written as follows: 

DK"^ 1 1 

— — = — 5^7^'fe'•Ji^^ + F" = -/^^7^'fc%/^^ + ^^^ (32) 

Dt 2m m 

T ^ <-) 

We by purpose use the symbols K"^ , K"^^, not P", S"'', because the latter sym- 
bols are reserved for canonical linear momentum and affine spin in contravariant 
representation. In the very special case of Lagrangian L = T—V{x^, s^a)i where 
T is given as in ipO]) . it is really so that these concepts are mutually identified 
via Legendre transformation, i.e., 

pa ^ ^abp^ = X° = mV, Y."^ = Y^cg'^ = K"^- 

However, for more general Lagrangians and for non-d'Alembertian models of 
the kinetic energy T, the above relationships become false, whereas ([32]) . (I33p 
remain true. 

The power, i.e., the time rate of work, is given by the following formula 
obtained by analogy with the mechanics of extended bodies [371 SS] • 



int ■ 



We can consider internal affine dynamics subject to additional constraints, just 
as in the case of extended affine bodies. First of all, let us consider gyroscopic 
constraints, i.e., assume the moving frame e to be permanently (7-orthonormal. 
Then, as mentioned above, U, is permanently g-skew-symmetric (17 is perma- 
nently 77-((5-)skew-symmetric). 

We assume here the validity of the d'Alembert model of constrained dynam- 
ics. Therefore, our equations of motion remain valid when on their right-hand 
sides we introduce some extra reaction forces responsible for keeping the con- 
straints, but, roughly speaking, not influencing the along-constraints motion. 
Let us denote the corresponding expression by Ff^, . They are to be added to 
the "true" applied dynamical quantities iV"''. According to the d'Alembert 
principle they are completely passive controls, i.e., they do not do work along 
any virtual motion compatible with constraints, i.e., 

Vb^ = gabF^V'' + gbcN%'^\ = 

for any V and for any g-skew-symmetric Vt. This means that Ffj — and N^^ = 
N'^, i.e., Fji vanishes and Nji is symmetric. Therefore, the effective, reaction- 
free system of equations of motion consists of p2p and the skew-symmetric part 
of ([33]) with algebraically substituted gyroscopic constraints gijC^ a^'' b = Vab = 
Sab, i.e., 

DfC"^ 1 



Dt 2m 
DS"'' 
Dt 



Af"", (34) 
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where, obviously, 

is the internal angular momentum (spin) and 

is the skew-symmetric moment, torque, i.e., generalized force coupled to rota- 
tional degrees of freedom. In this way one obtains the system of n{n + l)/2 
equations of motion imposed on the time-dependence of n{n -f l)/2 degrees of 
freedom of the infinitesimal gyroscope moving in M {n translational degrees of 
freedom and n{n — l)/2 gyroscopic ones). 

Returning to the explicit description in terms of the configuration variables 
we obtain that 

r)2 b rp-p°-, 
^^5i^ " Df^ " =Ar -TV . (35) 

In these equations only given (reaction-free) forces and moments are present 
and obviously 

In analogy to extended affine bodies one can also consider internal dynamics 
with other kinds of constraints. The most natural ones are those with the lucid 
group-theoretical structure. For example, for incompressible bodies, (|3T|) . ([33]) 
are to be replaced by their g-traceless parts: 

a ^ <i L jKL ^ „ ^ ^ L ^ab jKL _ ^jab „ Ajcd^ab 

i.e., 

R (x-b _ lg^^j^cd^ab\ ^ ^ab _ ig^^Ncdgab 

Dt \ n J n 

De-_KD^ _ 1 De-j^De^ ab 
Dt Dt Dt Dt ^ ■ 

Incompressibility constraints may be described analytically by equations 

det [e^i] = ^Jdetlffi^, 
or infinitesimally in any of two equivalent forms 

n\ = 0, h^A = 0. 
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Let us recall that gyroscopic constraints in Ricmannian space may be also de- 
scribed in terms of equivalent infinitesimal conditions: 



If only purely dilatational internal motion is admitted, i.e., 



or, equivalently. 



- -r!\.(5^ = 0, n^B - -^^c^^B = 0, 
n n 

then the internal motion is (on the basis of d'Alembert principle) completely 
described by the g-tiace of (|3T|) or ([33)) : 

gabe'^K^^^j"'' = 9abN^\ (36) 

9ab^^ dab^^^J +9abN . (37) 

If the body is subject to Weyl constraints, i.e., its internal degrees of freedom 
undergo only rigid rotations and dilatations, then equations of internal motion 
are given by ([M)) and ([57)) . or equivalently (pSjl and 

A very interesting example of constraints is that of rotation-less motion in 

M, 

n',=n/^g,kg''n^. (38) 

Then the d'Alembert principle implies that equations of internal motion (effec- 
tive ones, free of reaction forces) are given by the symmetric part of pip or 
dSSl), i.e., 

a ^ ^ L jKL I h ^ L jKL _ Arab , j^rba 

i.e., 

— (K"'' + X''") = -Pe°g De^L kl , De\ De^L kl , j^ab , nrba 
Dt^ ' Dt Dt Dt Dt -r -r . 

We have mentioned that when the geometry of AI is non-Euclidean (curved), 
then in a natural way aholonomic velocities and other aholonomic concepts 
appear. In the last example situation is even much more complicated, because, 
as mentioned in our papers about extended affine bodies [36l |37] , constraints 
(j38p are very essentially aholonomic even in mechanics of extended affine bodies 
in flat spaces; the more so in a non-Euclidean manifold. 

From now on we concentrate on the mechanics of infinitesimal gyroscopes 
and infinitesimal affine bodies without additional constraints, but very often 
with the special and strong stress on the mutual coupling between rotational 
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and deformative motion. This induces us to develop certain analytical proce- 
dures. They are even more than analytical procedures, because the underlying 
geometric techniques are interesting in themselves and simultaneously they re- 
veal certain very interesting mechanical facts. 

So, let us begin with the special case of infinitesimal gyroscope in a Riemann 
space (M, {}, g), where equations of motion are given by ([M)) or, more explicitly, 
by 

Let us stress an important fact. In mechanics of unconstrained infinitesimal 
affine bodies, equations of motion ([HO]) . (PTjl or ([55]) are directly apphcable 
even in the purely technical sense, because (a;*,e*A) are "good" independent 
(unconstrained) generalized coordinates. It is also so in expressions (^0]) . ((7f)) 
for the kinetic energy. However, after the gyroscopic constraints are imposed, 

9ije\e^ B VAB Sab) , 

the quantities b'a are no longer independent and cannot be used as general- 
ized coordinates. Even when we study the mutual coupling of rotational and 
deformative degrees of freedom, the quantities e^A are inconvenient, although 
well defined as generalized coordinates. In flat spaces the procedure is obvious: 
the system [c^a] as a matrix is subject to various polar, two-polar |49| . and 
similar decompositions, and then the mutual interplay between rotations and 
deformations is easily treatable. Something similar must be done here, but the 
direct methods based on the flat-space geometry are not applicable any longer. 
Instead, some geometric techniques based on orthonormal aholonomic reference 
frames may be developed. In the special case of gyroscopic constraints one 
can use then various well-known coordinates on the special orthogonal group 
SO(n,M) as a subsystem of well-defined generalized coordinates on {FM,g). 
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So, let E be some preestablished, fixed once for all fields of linear orthonormal 
frames on M, i.e., a cross-section of the subbundle {FM, g) C FM over M . By 
the way, the orthonormality demand is a technically simplifying one, but in 
certain problems some general frames would do. But let us assume E to be 
orthonormal in the (/-sense, i.e., 

g{EA, Eb) = gi-jE\E^B = Vab = Sab- 

The dual co-frame £'=(..., E"^, . . .) will be also used, 

{E"^, Eb) = E^iE^ B = 5^ B- 

It is obviously orthonormal with respect to the reciprocal (contravariant) metric 
Tj on M, 

E^.E^^^g'^ = ry^^ <5^^) . 

Remark: this is a field of frames, and it is defined all over in M and kept 
fixed. The natural question arises as to the choice of this aholonomic frame in 
M. In general the particular structure of the Riemann space {M,g) suggests 
some choices, both technically convenient and geometrically lucid, which are 
well suited to the problem. 

So, M is "inhabited" by continuum of orthonormal frames. The metric field 
on M may be expressed as follows: 

g = jfB^^ ^E"" = S^'^Ea ® E"". 

In the course of time the body moves in M. At the time instant i g M 
it is instantaneously placed at the geometric point x{t) € M and its internal 
configuration is then given by e{t) G F^^^^^M. Obviously, the vectors e^(t) € 
T3.(t)M, A = 1, n, may be expanded with respect to the just passed frame 
Ex{t) e f;(t)M, i.e., 

eA{t) = E^^^t^BR"" A{t). 
where, in the gyroscopic case, the matrix R is orthogonal, 

5ab = ScdR'^aR'^b, (39) 
or if for any reason we use a more general (non-standard) basis in R", then 

VAB = rjcoR^ aR^ B- 
Roughly, (|39p may be written in the usual matrix terms: 

R^R — In, 

where the n x n identity matrix is on the right-hand side. In this way the 
configuration space F{M,g), or more precisely F^{M,g) (orthonormal frames 
positively oriented with respect to some fixed standard of orientation) is iden- 
tified by the field E with the product manifold 

M X SO(n,R). 
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Obviously, this difFeomorphism depends explicitly on the choice of E. And of 
course everything is based on some topological assumption about M, namely, 
that it admits a globally defined smooth field of frames (and then, obviously, 
infinity of them). Instead of describing motion in terms of time-dependent quan- 
tities x^{t), e^A{t), we describe it in terms of quantities x^{t), R^B{t)- In other 
words, the motion is represented by a pair of curves: the one in M, i.e., the 
translational motion, and the other in SO(n, M)-internal rotational motion. And 
let us repeat that this trivialization (splitting) is possible only when the refer- 
ence field E is fixed and depends explicitly on it. One of the main advantages is 
a possibility of introducing good generalized coordinates instead of redundant 
quantities (x^e^^). Namely, SO(n,R) admits many well-investigated and ge- 
ometrically convenient coordinatizations like, first of all, canonical coordinates 
of the first kind (components of the rotation vector if n = 3), canonical coor- 
dinates of the second kind, and also some other possibilities. For example, in 
the physically interesting special case one uses also Euler angles or spherical 
variables in the space of rotation vectors. So, in any case we have at disposal 
various natural choices of "good", i.e., not redundant generalized coordinates 
(..., a;*, .. .), where fJ. = 1, n{n — l)/2, are the mentioned coordi- 
nates on SO(n, R). 

The same procedure may be applied in the general case of the infinitesimal 
afBnely-rigid body. The prescribed field of reference frames E = (• . . , Ea, ■ ■ ■) 
may be in principle quite general, but usually it is chosen as (7-orthonormal, even 
if deformations are admitted. The reason is that such a description is physically 
convenient when dealing with mutual interactions of rigid rotations and finite 
or infinitesimal homogeneous deformations. Just as previously we expand 

CAit) = E,f^t)B^''Ait), ^^B{t) = (E^, esit)) , 

but now ip{t) G GL(ri,M) is a general non-singular matrix, or usually a general 
positive-determinant matrix, ip{t) G GL+(n,M). In analytical procedures de- 
scribing realistic and simple isotropic problems we often represent (p in terms of 
the polar or two-polar decomposition [4^ . 

In this way the bundle of linear frames EM [E^M) is represented by a 
trivialization M x GL(n,R) {M x GL~''(ri,R)), also explicitly dependent on the 
choice of the frame E. Analytically, configurations are parameterized by gen- 
eralized coordinates (..., a;*, , (p^B^ ■ ■ ■)■ Let us stress once again that 
now it is not logically necessary, because the quantities (..., x\ , a, ■ ■ ■) 

themselves are good, independent generalized coordinates. The point is only 
that that, as mentioned, in many realistic problems generalized coordinates 
(..., a;*, , 95"^ B, . . .) are geometrically and technically more convenient and 
computationally effective. The more so it is when (p^ b are replaced by coordi- 
nates appearing in the polar and two-polar decompositions of ip. 

Let us begin with the general description of aflinely-rigid body in terms of 
prescribed reference frame E and later on consider the special case of infinites- 
imal gyroscopes. This will be also convenient for making the explicit use of 
the polar and two-polar decompositions, in particular, for discussing integrable 
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and superintegrable models in a two-dimensional space. The latter problem, 
i.e., two-dimensional infinitesimal gyroscopes or homogeneously deformable ones 
sliding over two-dimensional curved manifolds may be useful in the theory of 
microstructured (micropolar or micromorphic) plates. 

Obviously, the primary dynamical variables are the moving-frame vectors 
eA{t), or equivalently the moving- frame covectors, e'^(f), whereas the reference 
frames E are only the auxiliary quantities, does not matter how important for 
computational purposes. Let us express 

CAit) = EB{xit))ip^A{t), e^{t) = ip-'^B{t)E^{xit)), (40) 

where, obviously, e{t) e F^(^t^M, i.e., x{t) = TT{e{t)). 
We shall also need inverse formulas, i.e., 

EA{x{t)) = eBit)ip-'''Ait), E^{x{t)) = ^^B(t)e^(t). (41) 

To avoid unnecessary crowd of symbols, when it is not confusing, we shall omit 
the arguments t, x{t) in the above quantities. 

To express explicitly afhne velocities we must find formulas for the covariant 
derivatives DeA/E>t in terms of the quantities E, ip. 

The above expressions (|40l) imply that the covariant differentiation of M 9 
t 1-^ ca & EM along the curve of translational motion M. 3 t i-^ x{t) e M is 
given by 

DCA D , g . DEb B , IT '^'^^A 

—^ = -^\Eb^ a) = ^ a+Eb—, — . (42) 

Dt Dt^ ^ ' Dt ^ dt ^ ' 

Let us express the along-curve differentiation of Eb through its field-differen- 
tiation (well defined because E is defined as a field all over in M), 

DEb , ^ dx^ , s n dx^ , , 

^ = (V.i..)- = (Voi..)i.-.-, (43) 

where S/cEb is an abbreviation for VecE!b, and let us remind that for any 
vector field Y and tensor field T, VyT denotes the covariant derivative of T 
along Y. Analytically, 

VyT = Y'V^T. 

It is convenient to use aholonomic coefficients of our affine connection with 
respect to the field E, 

Vc'Eb = T^bc'Ea- (44) 

The usual holonomic coefficients of F with respect to coordinates x^ are given 
by the expression 

i jk — -ty A^ BCi^ ji^ k + J^A-Cy j,k, 

where the comma as usual denotes the partial differentiation. The second term 

r[EYjk E'AE^j,k (45) 
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denotes the teleparallelism connection induced by E. It is defined as follows: 

ViriE]]EA = 0, (46) 

i.e., it is the only afiine connection with respect to which all the fields Ea (of 
course, also their dual co-fields E"^) are parallel. Its curvature tensor vanishes 
and the torsion equals 

SlEYjk = E^AE^y^k^ = -j^E^A {E^i,k - E^k,j) ■ (47) 

The parallel transport with respect to r[£^] is path-independent. A tensor field 
T is rfSJ-parallel, V[p[£;]]T = 0, if its aholonomic components with respect to 
E are constant in M, 

T"'-"" h...ji = T^'-^" j,...j,Ei, (g, ■ ■ ■ (g, Ei^ (g, E'^' (g, ■ ■ ■ (g, E'^' , 
T^'-^'j^.-.J, = const. (48) 

If the field of frames E is holonomic, in particular, if it is given by the system 
of tangent vectors d/dx^ of a;*, then 

which means that we obtain the usual components of F with respect to local 

coordinates x*. 

Obviously, the quantity 

jk - J- [-c/j jk — -C' At BC^ k, 

i.e., in the coordinate- free form 

r - T[E] = T^bcEa <gE^<gE'=' 

is a tensor field, once contravariant and twice covariant, as it is always the case 
with the difference of affine connections. 

The teleparallelism torsion is directly related to what is well known in dif- 
ferential geometry as the aholonomic object of E, i.e., 

S[E] = ^n^BcEA ® ® E^. 

In other words the doubled _E-co- moving components of S[E] coincide with the 
Schouten aholonomic complex i7; the latter may be defined in terms of Lie 
brackets of basic vector fields, 

n^BC = 2S^BC = {E^, [Eb, Ec\) , (49) 

i.e., 

[Ea, Eb] = ^^abEc, dE^ = ^n^BcE^' A (50) 
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After some calculations performed on (gl]) with the use of (gOl), (SIl), (SSI), (Sll), 
we finally obtain the formula 



where the co-moving afHne velocity SI may be expressed as follows: 

oB T-F ,n ,„C f}E I ,-lB dtp'-'' A 

" A = Fi DC<y5 A"^ + (/? C , (51) 



where, as previously, 



at 



denotes the co-moving components of translational velocity. In this way the 
quantity 

^^^^ 



Dt / ^' ' Dt 

is expressed through the matrix ip and its time derivatives. This representation 
becomes very convenient when gyroscopic constraints are imposed. Let us ob- 
serve that in addition to the second term Lp~^dip/dt familiar from the flat-space 
theory, the above expression (|5ip contains the first term explicitly depending on 
the geometry of M (more precisely, on (M, F)) and on the choice of the auxiliary 
reference frame E. Roughly speaking, the quantity ip~^dif/dt is the co-moving 
affine velocity of internal degrees of freedom, as seen from the point of view of 
the just passed frame -Ea;(t)- And the first term represents the contribution to 
fl coming from the rotation-deformation of E itself. 

It is interesting to rewrite the equations (HH), dSO]), and §3^ 
in purely co- moving terms. After some calculations we obtain that 

dt 

+ ^h^DR''c^BV'' +F^, (52) 

= -n^Dn^cJ^'^' + N^''. (53) 

Obviously, the torsion here is taken into account, so strictly speaking they are 
CO- moving counterparts of ([28]) . (f29|) . 

Remark: do not confuse the three-index torsion tensor with the two-index spin 
one. The two-index quantity in the curvature term is vorticity. 

Using explicitly the balance form of the equations of motion for the linear 
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momentum and afRne spin we obtain that 



at m 



dK 



AB 



dt 



2m 

= -i?^^Jci3X^^+iV^^, (55) 



where, as previously, JbcJ^^ = 5^b- 

Obviously, we are allowed to use the usual time derivatives at the left-hand 
side of equations because P"*, K^^^ , V^, fi^^s are scalars from the point of view 
of geometry of M (although they are tensors in the micromaterial space R"). 

These are, so to speak, afiine Euler equations. By the way, when we impose 
metrical constraints, i.e., assume that permanently 

gieA,eB) = gije^e^ B = Sab, (/? G SO(n, M), (56) 

and according to the d'Alembert principle the reaction-free equations of internal 
motion, i.e., the skew-symmetric part of (|55p is only left, then we exactly obtain 
Euler equations for infinitesimal gyroscope in a non-Euclidean space. 

Let us observe that all the above equations of motion, in particular ([28]) . 
([29|) and the expression for the kinetic energy ((20|) . or its equivalent form ([2T|) . 
([2H) may be expressed in terms of (ISTI) . and in the case of gyroscopic motion 
this becomes technically unavoidable. The point is that a are not then inde- 
pendent generalized quantities, just due to the constraints (I56p . And although 
equations of motion in their general balance form (|34p are correct and free of 
reaction forces, in particular dynamical problems they are not directly useful 
in analytical calculations and analysis of the phase portraits. The quantities 
(fi^ L are not generalized coordinates either, because the matrix (p is orthogonal. 
However, as mentioned, ip € SO(n,R) may be easily parameterized on the basis 
of group-theoretical considerations, and the corresponding parameters are just 
proper generalized coordinates to be effectively used in equations of motion. Let 
g", a = l,n{n— l)/2, be such group parameters. Then taking some indepen- 
dent subsystem of ([M)) . e.g., one given by a < 6 (or conversely), we obtain a 
system of n{n — l)/2 second-order differential equations imposed on the time- 
dependence of n{n — l)/2 generalized coordinates. In variational problems it is 
more convenient to substitute the expressions (p{q) directly to the formula for the 
kinetic energy ([20|) or its Hamiltonian version ([25|) . and then obtain directly the 
corresponding Euler-Lagrange equations or (more convenient) Hamilton ones 
formulated in terms of Poisson brackets. Such an approach is better suited to 
the study for integrability problems and action-angle variables. Geometrically 
the most suggestive choice of generalized coordinates q" is that based on the ex- 
ponential representation of ip. The quantities are then canonical coordinates 
of the first kind on the group manifold of SO(n, R). More precisely, let Mkl & 
SO(n,R)' be basic skew-symmetric matrices; they have only two not vanishing 
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entries, namely ±1 in the rows and columns labelled by K,L [55^. Then R G 
SO(n,R) may be represented as follows: 

R{e) = exp (^Ie^^Mlk^ , 

where e is an arbitrary real skew-symmetric matrix, e^^ = — e^^, and, e.g., 
e^^ with K < L may be chosen as independent coordinates. The corresponding 
Mkl, K < L, form a basis of the Lie algebra SO(n,M)'. 

In the special case of n = 3, due to the exceptional isomorphism between 
axial vectors and skew-symmetric second-order tensors in M"^, it is customary to 
write [49] 

R (fc) = exp (fc^M^) , 
where Ma, A = 1, 2, 3, are basic skew-symmetric matrices, 

(e is the totally skew-symmetric Levi-Civita-Ricci symbol normalized by £123 = 
1), and are components of the so-called rotation vector in R'^. Its length 

k = ^ {k^f + {k^f + (k^f 

satisfies fc < tt if all possible direction versors n = k/k are admitted. Obviously, 
for any versor n we have 

R (nn) = R (— 7rn) . 

Therefore, the group manifold is represented in the space of M'^-vectors k by the 
ball of radius tt, with the proviso that antipodal points on the limiting sphere 
k — TT are identified. This provides a very nice form of understanding that the 
group S0(3,R) is doubly connected. 

In the three-dimensional case one uses also other parameterizations, depend- 
ing on dynamical details of the considered problem. For example, in certain dy- 
namical models of the spherical rigid body it is convenient to use the spherical 
variables (fc, (j>, 0) in the space of the rotation vectors fc [49] : 

k^ — k sin 9 cos 0, k^ ^ k sin 9 sin 0, k^ = k cos 9. 

In the theory of a free or heavy symmetric top (/i — I2, but in general 
Ii h), it is convenient to use Euler angles {(p, z9, V'), 

19,-0] = exp((^Af3)exp(z9Mi)exp(V'M3). 

One can also use canonical coordinates of the second kind {a,(3,"f), i.e., 

R[a,P,j] = exp(aMi)exp(/3M2)exp(7M3). 

Surprisingly enough, for the spherical rigid body (/i = I2 = ^3) the kinetic en- 
ergy expressions have almost the same form in coordinates (ip, i9, ip) and (a, 7). 
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If n = 2, the situation is much more simple, because then S0(2,R) is one- 
dimensional (thus, obviously, commutative) and R{ip) = exp((^e), where en = 
£22 = 0, ei2 — —621 = 1, i.e., explicitly 



In the case of rigid motion the second term in (|5ip . ip^^d(p/dt is always skew- 
symmetric, i.e., an element of the Lie algebra SO(n,M)'. More precisely, if we 
use artificial not orthogonal coordinates in R^, it is 77-skew-symmetric, i.e., 



saw, the total Q is also skew-symmetric (7/-skew-symmetric) and therefore, so is 
the first term of ([5T|) . By the way, the aholonomic connection coefficients r-'^^ic 
are then also skew-symmetric (more precisely, 77-skew-symmetric). Just as in the 
general affine motion, we can say that O is then the co-moving angular velocity, 
(p~^d(p/dt is the angular velocity with respect to the pre-fixed aholonomic frame 
E, and the first term is, in a sense, the angular velocity with which E itself 
rotates along the trajectory of the structured material point. 

In both the general affine and constrained gyroscopic motion the expression 
of kinetic energy through the above quantities becomes relatively lucid and com- 
putationally effective when the manifold (M,T,g) has some special structure, 
e.g., if it is a constant-curvature space and if the auxiliary field of frames E 
is appropriately chosen. The proper choice depends on the particular geome- 
try of M, and it is also a matter of some inventive intuition. In such special 
cases the problem may be effectively studied on the rigorous analytical level 
and, in particular, interesting results concerning integrability and degeneracy 
(superintegrability) may be obtained [Tl[Tl[Iil[Tl[ni|ISl[Mlll51[2Sl[lI]. 

Unless otherwise stated, from now on we concentrate on Riemannian spaces, 
when r is the Levi-Civita connection built of g. The general case, when T,g 
are unrelated may be also interesting in itself, but certainly there are some 
difficulties to be overcome, because even in the rigid motion the afiine velocity 
fl is not g-skew-symmetric and its co-moving representation 17 is not r]-{S-)skew- 
symmetric. The fact that they do not belong to the corresponding Lie algebras 
of the rotational groups SO{TxM, g^), SO{n,M.) obscures their interpretation. 
If (M, r,g) is a Riemann-Cartan space but F is not symmetric (i.e., it is not 
Levi-Civita connection), then affine velocities are skew-symmetric (respectively 
in a g- or ?7-sense), but they are explicitly dependent on the torsion tensor. This 
also leads to certain interpretation problems. Namely, one obtains different 
expressions for the kinetic energy and different equations of motion when one 



i?(^) 



cos (p — sm ip 
sin ip cos (fi 




If the connection T is metrical, Vpjff — (Riemann-Cartan space), then, as we 
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in the definition of angular velocity. 

In the general case of affine motion the formula (|51[) will be written in any 
of the abbreviated forms 



where 



n,,^B^n{ri)\ - ^"'^c:^. (58) 

The labels "dr" and "rl" refer respectively to "drift" (or "drive") and "relative". 
The reason is that, as mentioned, ip refers to affine rotations with respect to the 
just passed prescribed reference frame E; the first term describes the time rate 
of affine rotations contained in the field E itself. When gyroscopic constraints 
are imposed, all these fi-objects become skew-symmetric angular velocities. To 
stress this, sometimes, but not always, we shall use then the symbols cD, cDdr = 
a)(dr) and lu^\ — u)(t\). 

One of analytical advantages following from the prescribed reference frame 
E is the possibility of using the polar and two-polar decompositions [49] , 

ip^UA^BU ^ LDR-\ 

where ?7, L, R are orthogonal (more precisely, r^-orthogonal when artificial non- 
orthonormal coordinates are used in R"), A, B are symmetric (//-symmetric), 
D is diagonal, and obviously, 

B = UAU-\ 

As usual, U, L, R denote fictitious gyroscopic degrees of freedom extracted 
from ip £ GL(n,R) j49]. The corresponding "co-moving" angular velocities are 
given by the expressions 

^ 1 dU ^ ^ T dL ^ „ 1 dR , , 

-r.^U-^-, X. = L-^, ^. = i?"V. (59) 

Obviously, as usual, the "spatial" representation may be used: 

fi?7 1 dL , dR , 

However, in calculations appearing in practical problems the "co-moving" ob- 
jects are more convenient. Obviously, in the two-dimensional world, when n = 2, 
these representations coincide. 

After some calculations one can show that the kinetic energy of internal 
motion Tint (|20p may be expressed in the following way in terms of the polar 
decomposition: 

Ti„t = -iTr(AJAS2)+Tr('^J^C)+^Tr( J(^) ), (60) 
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where 

a) = cDdr + iirf = Sdr + t^"^-^, (61) 



and obviously ujdi is the restriction of fidr (|57|) to the ?7-rigid motion: 

Sd/s = [/"i-^f r^zjcf/'^Bt/^BF^. (62) 

The formula (j60p is written in the standard orthonormal coordinates in M". 
Otherwise, when tiab is admitted to be different than Sab, we have to replace 
the symbol J in ((60|) by J^, where 

^„^B J'^^VCB- 

Geometrically, J,, is a mixed tensor in R", i.e., linear endomorphism of M", G 
L(n,R), whereas J itself is a twice contravariant tensor. 

The two-polar decomposition becomes analytically useful in doubly-isotropic 
dynamical problems, i.e., ones isotropic both in the physical space M and the 
micromaterial space. This double isotropy imposes certain restrictions both on 
the kinetic and potential energies. What concerns the very kinetic energy, the 
inertial tensor must be proportional to 77, i.e., 

jAB T AB J A TxA 

Then one can show that (l60l) becomes 



where now 



Tint = - ^Tr (D^x^) - ^Tr (z?^^) + /Tr (^DxDd^ , (63) 



5 = «-f . (64) 

X = Xdr + Xrl = Xdr + i , (65) 

Xdr^s = L-'^fT^dcL''bL''eV''. (66) 

The last formula is quite analogous to ((62)) . Just like there, x contains the 
"drive" term built of the connection coefficients. It is only the L-rotation that is 
coupled in this way to spatial geometry; the i?-rotation is geometry- independent. 

Again we conclude that (|63p is structurally identical with the corresponding 
formula for extended affine bodies [37] with the proviso however that x contains 
the drive-term. The expression for 1} is free of such a correction. Everything 
that has to do with (M, F, (ji)-geometry is absorbed by the x-term. 



5 Special two-dimensional problems 

We shall consider now some special two-dimensional cases, i.e., when n — 2. 
Therefore, for the infinitesimal rigid body (infinitesimal gyroscope) we are deal- 
ing with three degrees of freedom: two translational and one internal, rotational. 
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If no gyroscopic constraints arc imposed and the internal motion is affine, then 
of course there are four internal degrees of freedom; together with translational 
motion one obtains six degrees of freedom. The resulting models are interesting 
in themselves from the point of view of pure analytical mechanics, in partic- 
ular, some integrability and hyperintegrability (degeneracy) problems may be 
effectively studied. Obviously, the explicit analytical results exist only in Rie- 
mann manifolds {M, g) with some peculiar structure, first of all (but not only) 
in constant-curvature spaces. Some practical applications of two-dimensional 
models also seem to be possible, e.g., in geophysical problems, in mechanics of 
structured micropolar and micromorphic shells, etc. What concerns geophysics, 
we mean, e.g., motion of continental plates. Motion of pollutions like oil spots 
on the oceanic surface is another suggestive example. 

Below we consider in some details three kinds of two-dimensional problems, 
namely motion of structured material points on the sphere, pseudo-sphere (Lo- 
batchevski space), and torus manifolds with geometry induced by injections in 
the three-dimensional Euclidean space M^. Many interesting dynamical models, 
including some quite realistic ones, may be effectively investigated in analytical 
terms. It is not very surprising in spherical and pseudo-spherical geometries 
because of exceptional properties of constant-curvature spaces. But there exist 
also nice completely integrable models on the M^-injected torus. Perhaps this 
may have something to do with that that all these manifolds are algebraic ones 
(of the second degree in the spherical and pseudo-spherical cases, and of the 
fourth degree in toroidal geometry). 

5.1 Spherical case 

Let us begin with the spherical geometry. The two-dimensional "world" will be 
realized as the two-dimensional sphere in M.^ with the radius R and the centre 
at the beginning of coordinates. Of course, the radius has the intrinsic sense, 
namely the scalar curvature equals 

We introduce the "polar" coordinates (r, tp) . Obviously, r is the geodetic 
distance measured from the "North Pole", i.e., § = r/R is a. modified "ge- 
ographic latitude", and (p is "geographic longitude", i.e.. Rep is the distance 
measured along "equator". So, ip, are usual polar angles in M^, and R is 
the fixed radius-distance from the origin. Euclidean coordinates of points on 
S^{0, R) e are given by the usual expressions: 

X = Rsin'&cosip, y = Rsm'&sinip, z = Rcosi}, + + = R? . 

Obviously, r runs over the range [0, ttR] from the "North Pole" to the "South 
Pole" , and has the usual range [0, 27r] of the polar angle. Coordinates are 
singular at r = 0, r = -kR, and there is the obvious 0, 2-k ambiguity of the 
longitude p. 

Restricting the metric tensor of 9? (roughly speaking, the metric element 
dx'^+dy'^+dz'^) to S^{0, R) we obtain the obvious expression for the arc element 
on the sphere 

ds^ = dr^ + R^ sin'' ^dip^. (67) 
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So, the kinetic energy of translational motion is given by 

where m denotes the mass of the material point. 

Even for the purely translational motion some interesting questions arise, 
e.g., what are spherically symmetric potentials V{r) for which all orbits are 
closed? Obviously we mean problems based on Lagrangians 

itr =Tt,-l/(r). 

This is a counterpart of the famous Bertrand problem in R-^. And it may be 
shown that the answer is similar [5D1 [3S1 HO] , i.e., the possible potentials are as 
follows: 

(i) oscillatory potentials: 

l^(r) = |i?V^, (69) 
(ii) Kepler-Coulomb potentials: 

nO = -|ctg^- (70) 

Obviously, with the spherical topology also the geodetic problem belongs here: 

{in) V{r) = 0, i.e., (in a sense) the special case of {i) or (ii) when x — 0, 
a = 0. 

There is an obvious correspondence with the flat-space Bertrand problem; it is 
suggested by the very asymptotics for r « 0, i.e., 

V(r) « -r\ V{r) « --. 

Obviously, this is a rough argument, but it may be shown [35l |40] that there 
exists a rigorous isomorphism based on the projective geometry. 

The mentioned Bertrand models lead to completely integrable and maxi- 
mally degenerate (hyperintegrable) problems. But even for the simplest, i.e., 
geodetic, models with the internal degrees of freedom the situation drastically 
changes. There exist interesting and practically applicable integrable models, 
but as a rule interaction with internal degrees of freedom reduces or completely 
removes degeneracy. 

Let us begin with the gyroscopic model of internal motion. Unlike the gen- 
eral case, in two-dimensional problems with the constant-curvature spaces, gy- 
roscopic problem is simpler than affine one. Moreover, it very simplifies study 
of the affine case. 
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The first step is to introduce an appropriate field E, i.e., the auxiliary and 
fixed once for all orthonormal aholonomic reference frame. This is often the 
matter of inventive guessing. In this case it is natural to expect that the natural 
base (holonomic one) tangent to coordinate lines is a good starting point. It 
consists of the vector fields 



d_ 

dr' 



d_ 
d(p' 



Obviously, we mean here the well-known identification between vector fields and 
first-order differential operators, i.e.. 



X = X' 



d 



(X^ are components with respect to the manifold coordinates x^). 

Using more traditional language: if coordinates are ordered as (r, 99), then £r 
and Sip have respectively components [1, 0]"^ and [0. 1]"^. This holonomic system 
is not appropriate however, because it is not orthonormal (no holonomic system 
may be so in a curved Riemann space). However, it is evidently orthogonal (r, 
(p are orthogonal coordinates): 



9r(p — 9cpr — O5 



1, 



R 



The lengths of f^, £ip are obviously given by \\£r\ 
it is natural to expect that the normalized fields 



1, \\£. 



Rsm{r/R). So, 



or 



f 



d 



R sin dip 



(71) 



will form a convenient orthonormal frame E. Let us stress that E is aholonomic, 
although directions oiEr, E^ are tangent respectively to coordinate lines of r, ip. 
The reason is just the "crossed- variables" normalization, i.e., £^ is multiplied 
by a function dependent on r. There are no coordinates r', ip' for which we 
would have E^ = d/dr', E^ = d/dip'. The orthonormal frame e = (61,62) 
describing the internal configuration is obtained from the fixed aholonomic frame 
E = [Er, E^) with the help of some time-dependent orthogonal matrix ip = U, 
i.e., 



61 = E^U\+E2U^i=ErU\+E^U'^i, 

62 = EiU^2 + E2U^2 = ErU''2 + E^U'^2. 

Obviously, U may be parameterized in the usual way, i.e.. 



U = 



cos "0 — sin -0 
sin tp cos if) 



and 





■ -1 ■ 


_ 


■ -1 " 


(Drl = Wrl 


1 


~ Itt 


1 
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Remark: do not confuse the scalar quantity ujr\ with the second order skew- 
symmetric tensor uj^-i parameterized by it. 
Similarly we write that 



-1 

1 



-1 

1 



where 



cos ■ 



r dip 



cos-d 

dip 
— 4 
dt 



dip 
It' 



r dip 



dijj 

H 



COS'S 



dip 
It' 



(72) 



Obviously, the above expressions are the special cases of ([59)) . ((6T|) . ((62|) and 
follow easily from the formulas for Christoffel symbols (Levi-Civita connection) 
on the sphere, i.e.. 



R . r 
— sin — , 
2 i?' 



1 r 

r'^ — r''' — — cttr— 



rip ^ ipr 

(the remaining coefficients vanish). In the aholonomic representation: 

1 r „,„ 1 



r 



ipip — r'^^^ R' ^'^rip — Jl^^^ Jl' 



but let us notice that 

ipr ^ / ^ rip , 

i.e., the aholonomic coefficients are in general not symmetric. The remaining 
aholonomic coefficients also vanish. Obviously, the above V^bc are exactly the 
formerly used T'^bc] simply in some particular concrete cases this way of writing 
is more suggestive and intuitive. 

After simple calculations based on the formula (|20p or rather on its restriction 
to the gyroscopic motion, one obtains the expected expression: 



T 



i.e., briefly, 




(73) 



i^^sin^^- 
R 




T\-nt — 



2 ' 



where lu is given by (j72p . The scalar quantity /, i.e., the inertial moments of 
the plane rotator is related to the tensor J by the formula 



VabJ 



AB 



TrJ„ 



i.e., in Euclidean material coordinates (t]ab — 6ab) simply / = TrJ. In the 
absence of deformations the internal inertia is controlled only by this single 
scalar. This is the peculiarity of the "two-dimensional world" . 
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For certain reasons it will be convenient to rewrite the formula ([75]) in terms 
of the variable -d — r/R, i.e., 



T 




(74) 



— + sin^ i9 




2 \ dt dt 



It is seen that if formally {ip, d, ip) are interpreted as Euler angles (respec- 
tively the precession, nutation, and rotation), the above expression is formally 
identical with the kinetic energy of the three-dimensional symmetric rigid body 
(without translations) with the main moments of inertia given respectively by 



h=h= raW 



h = I. 



If / = mR^ one obtain the expression for the spherical top. 

There is nothing surprising in the mentioned isomorphism because the quo- 
tient manifold SO(3,M)/SO(2,M) may be identified with 8^(0, 1) (or with any 
S^(0, i?)) in a natural way. Projecting the motion of the three-dimensional 
symmetric top onto the quotient sphere-manifold we obtain two-dimensional 
translational motion; the one-dimensional subgroup of rotations about the z- 
axis refers to the internal motion of the two-dimensional rotator. 

The projection procedure is exactly compatible with the mentioned corre- 
spondence between Euler angles in S0(3,R) and our generalized coordinates 
d = r/R, ip of the infinitesimal rotator in S^(0, R). 

Let {/{(fijiifjip) G S0(3,i?) be just the element labelled by the Euler angles 
</?, ip, thus 

Uiip, ^, V) = U,{ip)U,{m.W, (75) 

where Uz, Ux are rotations respectively around the z- and x-axes; angles of rota- 
tions are indicated as arguments. Calculating the "co-moving angular velocity" 



K = U — — 

dt 

of this fictitious three-dimensional top one obtains that 



where 




0-1 
1 



Hi 
K2 
3<3 



X2 



1 


-10 



^3 



-1 





• ^ ■ I d-d 

smt/smw h cos w — , 

dt dt 

■ q i^f . ,d^ 

smv cosip— sm ■(/;—, 

dt dt 



cos t9 



dip dip 
H ^ ~dt' 



(76) 

(77) 

(78) 
(79) 
(80) 
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In expression we easily recognize ([7^ . i.e., the expression for the one- 
component angular velocity of the two-dimensional rotator. Calculating for- 
mally the kinetic energy of the three-dimensional symmetric SO(3,R)-top, i.e., 



(81) 



and substituting K ~ mR^, d — r/R, we obtain exactly ([73|) . i.e., ([74]). 

The principal fibre bundle of orthonormal frames F(S^,g) over the two- 
dimensional sphere (with its induced metric g) has S0(2,R) as the structural 
group and it may be itself identified with S0(3,M). This is just the essence 
of the above identification between S0(3,M) and its quotient manifold ~ 
SO (3, M)/S0(2, M). Some important invariance problems appear there. Namely, 
as usual in analytical mechanics, the kinetic energy (|73p . (|74p may be identified 
with some Rieniannian structure on the configuration space. Let us write down 
our kinetic energy in the following form with the explicitly separated mass factor: 



Just as above, our generalized coordinates g', i — 1,2,3, are the variables 
{r,ip,tp) written just in this direction. On the level of S0(3,M) as identified 
with F(S^,(7), they are equivalent to the Euler angles {'d,ip,ip), where = r/R. 
After some calculations one obtains that 











i?2 sin^ -g, + -L cos^ ^ -^-cos^ 

H m H m H 







In the special case / = mR^ one obtains that G simplifies to G, where 











i?2 i?2 r_ 

R^cos^ i?2 



(82) 



The corresponding expressions for the weight-one volume densities are as follows: 



, / r 
G = R\l — am—, 
m R 



'G = R^ sin ■ 



R 



The contravariant inverse metric G*-' {G^'^Gkj = S^j) is given as follows: 

1 







fi2 f 



COS 



i?2 



R mil „fo.2 r 
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and, obviously, 



1 















1 








fl2 sin^ ^ 



If we once identify the bundle manifold F{S^,g) with the group S0(3,R), then 
we can consider the action of two transformation groups. They are, obviously, 
represented by the left and right regular actions of S0(3,R) on itself, i.e.. 



where X,V & S0(3,R). From the point of view of the original configuration 
space F{S'^,g), these groups act both on the positions of material point in 5^ 
and orientations of the attached bases, i.e., internal configurations. 

The above metrics G are invariant under the total left-acting group SO (3, R). 
What concerns the right actions, in general they are invariant only under the 
subgroup S0(2,R) C S0(3,R) interpreted as the group of rotations about the 
z-axis in R'^. Therefore, G has the four-dimensional isometry group S0(3,M) x 
S0(2,R). Only in the special case / — mR^ the metric tensor G is invariant 
under S0(3,R) x S0(3,R), i.e., under all left and right regular translations. 
Therefore, on the level of S0(3, R)-description we are dealing then with the 
spherical top and up to a constant multiplier G becomes the Killing metric on 
S0(3,R). Of course, this special case is not interesting from the point of view 
of our primary model for an infinitesimal rigid body moving in 5^(0, R) because 
it is there physically too exotic and too exceptional. 

For the potential systems with Lagrangians L = T — V{q) the Legendre 
transformation has the usual form: 



X <-* XV, 



Inverting it. 




(83) 



one obtains the Hamiltonian H, i.e.. 




T=—G^^{q)p.p,. 
2m 
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As a rule, the presence of forces, i.e., not constant V{q), reduces strongly the 
aforementioned symmetry of geodetic models. 

We are interested here in integrability and hyperintegrability (degeneracy) 
problems, thus, we concentrate our attention on the Hamilton- Jacobi equation, 
i.e., 

its separability, and the action-angle variables. We are dealing only with time- 
independent problems, thus, S{t, q) is always sought in the following form: 

S{t,q)^-Et + So{q), 

where, obviously, the reduced action 5*0 satisfies the time-independent Hamilton- 
Jacobi equation, i.e., 

h{<i\^4^]=E. (84) 



dqi 

Let us write down explicitly a few formulas describing Legendre transforma- 
tion and kinetic energy for the rigid body moving in two-dimensional spherical 
world. 

Denoting canonical momenta conjugate to coordinates (r, 7p) respectively 
by {pr,Pip,P^), we have that 

Pr = mf, 

r 

P0 = lip + I cos— ip. 

H 

The resulting geodetic Hamiltonian has the following form: 

pf pI ~ "^PvP^ i + ^^"^^ i + Ti)pl 

pj P% - 2p<pP^ cos ti+pI mR^ - / 2 
2mi?2 sin^ ^ + 2mm I 

It is seen again how this expression simplifies in the special case when / = mE? . 

Without gyroscopic degree of freedom, when T is reduced to Ttr given by 
the formula (l68l). T is reduced to 



^2 „2 



Pr , Pf 



2m 2mi?2 gin^ ^ ' 
Then there exists the class of separable potentials 

V{r,ip) = Vr{r)+ ^^^^^ 



i?2 sin^ ^ ' 



54 



With rotational degree of freedom the metric tensor G is not diagonal in natural 
coordinates (r, if, ip) and there is no direct analogue of above integrable models. 
Of course, in three-dimensional Riemann manifolds there are always orthogonal 
coordinates, but it is fairly not obvious whether there exist diagonalizing coor- 
dinates admitting a physically reasonable class of potentials treatable in terms 
of the separation-of-variables method. 

It is seen however that ip are cyclic variables in the kinetic energy term; 
this focuses our attention on the models where the potential energy also does not 
depend on -ip. These angles are then cyclic variables for the total Lagrangian 
L = T ~ V and the corresponding conjugate momenta Pip,Pip are constants of 
motion. The resulting models, including geodetic ones = 0), are completely 
integrable and one can analyze them by means of the separation-of-variables 
method. 

As usual when dealing with cyclic variables, we seek the reduced action 5*0 
in the following form: 



where £, s are integration constants, i.e., the dependence of S^{ip) and S^^ip) 
on their arguments is postulated as linear. Together we have three integration 
constants E,£, s, just as it should be in a complete integral for the system with 
three degrees of freedom. As a matter of fact, due to the assumed symmetry, 
the problem reduces to the one-dimensional one for Sr and substituting (|85p 
into (j84p . we obtain the ordinary differential equation 



Therefore, Pr — dSr/dr equals plus-minus (depending on the phase of motion) 
square root of the rigid-hand side of the expression above. Obviously, this is 
well defined only in the classically admissible region between the turning points. 
It is clear that 



So{r, ip, ip; E, £, s) = Sr{r; E) + £ip + sip, 



(85) 




(86) 



Pf = 




dSa _ dS^ 
dip dip 



d^p ~ d^ ~ 



OSq dS^ 



The corresponding action variables are given as follows: 
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The radial action variable 

Jr = (f> Prdr 



equals the doubled integral of the square-rooted right-hand side of between 
the turning points, i.e., between nulls of (j86p . Substituting there 



we obtain the expression 



2m {E - V{r)) - ' ^ , + —^dr. 



Let us remind that the above expression, i.e., the contour integral of the dif- 
ferential one- form prdr along the corresponding orbit in the two-dimensional 
phase space of the (r,pr)-variables, equals 



Jr ~ 2 / pr{r)dr, 



where r^i^,r^s^y^ denote respectively the left and right turning points of the 
r-motion. 

When some explicit form of V{r) is assumed and substituted to then 
in principle the integral may be calculated and one obtains the functional de- 
pendence of the quantity Jr on the integration constants E, Jip^J^, i.e., 

Jr — Jr ( ; dip^ dtp') ■ 

Once calculated, this expression may be in principle solved with respect to the 
energy E, i.e., 

E — Ti {Jr, J^, J-p) ■ 

In this way the energy is expressed in terms of action variables. Substituting 
this expression and (j87p into (j85p , we obtain the generating function of canonical 
transformation from the original phase-space coordinates (r, Lp,il}]pr,Pif,Pti)) to 
the action-angle variables (0^, O;^, 0^; Jr, Jip, J^), i.e., 

Sq (r, (fi, Jr,J^, Jtp) = Sr {r,Ti{Jr, J^, J-4,)) + ^d^ + ^d,p. 

Ztt Ztt 

The resulting angle quantities are given as follows: 

« « _ dSo 

Odr Od^ OJtP 

This "angles" are meant modulo 1 but perhaps more intuitive are angles taken 
modulo 2tt, i.e., 

Qr=2TTQr, Q^^2nQ^, Qp^2TTe,p. 
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Let us observe that usually 



The fundamental frequencies are given as follows: 

_ m _ dn _ dn 

dJr ' dJy, ' '^^ 9 J,/, ' 

and their circular counterparts are the corresponding 27r-multipliers, i.e., 

As usual, the action variables are constants of motion, i.e., 

dJr _ dH _ dJ^ _ dH _ dJ^ _ dH _ 

~dt~^r~ ' ~dr~'d%~ ' "dT ~ 56^ ~ 

And as usual, the angular variables perform uniform "librations" because 

de' m i,^, 
-d^ = djr''^'^^ = 

thus, _ 

@\t) = v\J)t + a\ e\t) = oj\j)t + 13\ 

and, obviously, the index i runs over the labels r, ip,ip. 

A very important point is hyperintegrability, i.e., degeneracy. Let us remind 
that a multiply periodic system with n degrees of freedom (an integrable one 
with an open set of bounded trajectories) is said to be fc-fold degenerate (or 
(n — fc)-fold periodic) when fundamental frequencies satisfy a system of k (but 
not more) independent equations, i.e., 

n«iz/*(j)=0, a = T7fc, 

where are integers, i.e., G Z (but nothing changes when we say that they 
are rationals, i.e., n"i € Q). By "independent equations" we mean obviously 
that Rank [n^i] = k. Integers rj"; are to be constant, independent of J, thus, 
the same for all possible frequencies. If the above equations hold exceptionally 
for certain exclusive values of J, one says about accidental degeneracy. 

Of course the n-dimensional manifolds Ji = const are tori. While { J^, J^} = 
0, i.e., the action variables are in Poisson-bracket involution, these tori are La- 
grangian submanifolds, i.e., the symplectic two- form vanishes when evaluated 
on any pair of their tangent vectors at any point, and they are maximal con- 
nected manifolds of this property. If the system is /e-fold degenerate ((n — fc)-fold 
periodic), then the closure of every trajectory is an (n — A;)-dimensional isotropic 
torus (the symplectic form vanishes when evaluated at any pair of tangent vec- 
tors attached at any point). The system orbits are dense in these tori. The 
mentioned (n — fc)-tori form a regular foliation (congruence) of partially dis- 
joint submanifolds on the Lagrangian tori = const, a = l,n. For any fc-fold 
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degenerate system Hamiltonian depends on the action variables Ji , . . . , J„ in a 
very peculiar way, namely, it is a function of some (n — k) linear combination of 
Ji, . . . , J„ with integer coefficients (and of no less number of such combinations). 
The extreme examples are as follows: 

(i) a k = 0, then the system is completely nondcgencratc, any orbit fills 
densely some torus Ja = const, a = 1, n, and H depends on all Ja-variables 
in an essential way, i.e., there is no possibility to superpose them into a 
smaller numbers of combinations with integer coefficients, 

(ii) if fc = n — 1, then all bounded trajectories arc periodic, i.e., they are 
topological circles (one-dimensional tori) regularly foliating the tori = 
const, a = l,n, so that the quotient manifolds are (n— 1) -dimensional tori. 
Then there exists a combination of Ja-variables with integer coefficients, 
i.e., J = n" Ja, n" G Z, such that H depends on all Jo-s only through J. 

The first special case is completely not resonant and has n independent fre- 
quencies (periods). The second, quite opposite, case is maximally resonant, i.e., 
there is only one fundamental frequency (period) and trajectories are known 
from the elementary mechanics as Lissajous figures in the original configuration 
space. 

If the system is fc-fold degenerate, then the original action-angle variables 
obtained from the space-time coordinates q^,Pi, i = l,n, may be replaced by 
some new ones, better expressing the degeneracy (hyperintegrability) structure. 
Namely, it is always possible to introduce some quantities 9", Jo such that 

(i) Ja = JbN\, Q- = N\Q\ 

(ii) N"i) G Z, i.e., entries of the matrix N are integers, 

(m) det[iV%] = ±1, 

(iv) H = H ^Ji, . . . , Jn-k^ , i.e., Hamiltonian depends only on the first (n—k)- 
tuple of the action variables, 

(v) the first (n — A;) basic frequencies u^, r = 1, (n — k), are independent over 
integers, i.e., the equations 

^ ^ dJr 

r—l r—1 y-'^r 

imply that the (n — A:) x (n — fc) matrix rv' ^ vanishes, i.e., Z 9 n'*,. = 0. 

Obviously, the remaining k frequencies, i.e., = dTi/dJp = with p = 
{n — k + l),n, vanish and the corresponding new action variables 0^, p = 
{n — k + l),n, arc constants of motion. Together we have (n -|- k) constants 
of motion, global and smooth ones: 

I'l) • • • ) 'Jm , . . . , w . 
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What concerns globality, according to the above conditions (i) and (m), the 
quantities 8°, a = l,n, are "good" angular variables (they would not be so 
without these conditions). They are of course multivalued in a "harmless" way 
when property treated, just like the angular variables. To avoid this multival- 
uedness, one can replace them by unimodular complex numbers, i.e., 

C = exp {ie") . 

Let us go back from these general digressions to our special models, beginning 
from thus, to n = 3. 

It happens often, e.g., in the case of symmetry in the material point mechan- 
ics, that some partial features of the problem may be seen from the formulas 
like (j88p for even without calculating the integral and without even assum- 
ing anything about the shape of V{r). For example, in central problems of 
the material point motion, when the spherical coordinates r, ip are used, it 
is seen that the corresponding action variables , J^p enter the under-square- 
root expression through the combination + Jip, so the one- fold degeneracy is 
obvious from the very beginning. The particular shapes of V{r) are necessary 
for answering the question concerning total degeneracy. And then, according to 
what is known from somewhere else, really one can prove after the explicit in- 
tegration, that the total degeneracy occurs for the Kepler- Coulumb (attractive) 
problem and the isotropic harmonic oscillator, i.e., 

V = --, V^-r^, a > 0, >c>0. 
r 2 



In the formula ((88)) nothing like this is seen, even for the very special, highly- 
symmetric case / — mr'^ , i.e.. 



Jm — '^JtoJjb cos -t- Jl 
Jr =fxl2miE- V{r)) - ^ 4,.^,,,/. ' ^r. (89) 

Without performing some calculations, one does not see anything even in the 
geodetic case, when V = Q. But everybody knows that the problem is then 
completely degenerate because it is isomorphic, as mentioned, with the free 
spherical top, each orbits of which are closed (one-parameter group and their 
cosets). But without the explicit calculation nothing may be decided because 
Jip and do not combine integer- wise under the square root. Moreover, there 
are some unexpected strange technical problems when calculating then the re- 
lationship between E and J-variables. This is probably the reason one can 
hardly find the action-angle analysis in the rigid body mechanics. Passing to 
the isomorphic problem for the rigid body we obtain that 



, , J2 - 2 J^Jw, cos?? -t- Jl 

where "d = r/R. Substituting here x = — cost?, we obtain for the geodetic case 
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that 



j^ = -(f>d-2IEx^-^x + 2IE-:^^-^- ^ -. (91) 

^ " 27r2 47r2 (a;-l)(a; + l) ^ ' 

This integral may be calculated according to the usual rule for integrals of the 
form 

TZ ^\/ ax"^ +bx + c, x^ , 

where TZ denotes a rational function of indicated expression. Of course, the 
corresponding indefinite integral is elementary one and can be calculated, then 
the resulting definite integral (the doubled value between the turning points) 
may be in principle obtained in this way. There is however plenty of possibilities 
of making mistakes and it is much more convenient to use here the method of 
complex integration, as elaborated by Max Born in analytical mechanics and the 
old quantum theory. Replacing x by the complex variable z, we see that there are 
exactly two branch points, i.e., just the classical turning points. Therefore, the 
action variable (J,.. J^) may be calculated as the contour integral along the 
path infinitesimally surrounding the cut joining the branch point along the real 
axis. There are also three poles, i.e., z = —1, z = +1, 2: = oo. Surrounding these 
poles by additional contours and taking the composed unconnected contour 
consisting of the above four ones, we finally conclude that 

= — 27riRes_i — 27riResi — 27riReS(x), 

where, obviously, residua are calculated for the total integrand. After some 
calculations one obtains that 

Res_i = JjLpJlMi^ 

\Jif, + J1l}\ . 

Resi = — —I, 

47r 

Resoo = V2TEi. 

Thus, finally, 



47r\/2lB = 2J.a + \J^ - J^\ + | + J^-l • 

Only now, after all calculations, the degeneracy is seen on the level of action 
variables, although it was obvious from the very beginning from the spherical 

top analogy. There are however some delicate points, namely, the structure of 
this total degeneracy is a little bit different in the four regions of the phase 
space: 



{i) {Jip > J^) A (J<^ > — J^j,), i.e., > I J^l, then opening the absolute value 
sign we obtain that 

„ (Jtf + Jy) 
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(m) (J<p < J^) A [J^ > -Jtp), i.e., \ J^\ < J^, then 
(iii) (Ji^ > J^) A (J^ < -J^,), i.e., | < | J,/,] = -J^, then 

(if) (J<p < J^) A (J;^ < -J^,), i.e., Jip < - |J^|, then 

( - J^)'^ 
87r2/ ■ 



It is seen that there are four open regions of the phase space in which the formula 
fooks shghtly different. Nevertheless, the difference is not very essential. In any 
of these region one can introduce such action variables that only one of them is 
essential for the energy, i.e.. 

On the level of the old quantum theory according to the Bohr-Sommerfeld pos- 
tulates we have that 

and we obtain that 



J — nh, n ^TL^ 
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i.e., just as it should be. 

By the way, it is seen from the formulas O, (El) that the 

above class of integrable problems contains also interesting non-geodetic models. 
Obviously, this is the case when the structure of V is appropriately sited to other 
(geodetic) terms occurring under the square root sign. In any case, it is seen 
from (pij) that potentials of the type 

V{x) — ax^ + f3x, a, (3 = const 

do not violate the explicit solvability in terms of elementary integrals. If /3 7^ 0, 
then they also do not violate the total degeneracy because the term with the 
constant a will be simply absorbed by —2IE and will regauge it. In other words, 
in ((90| we can simply use that 



Vi^) =5003^19-^/3 cos I?, 

i.e., 

V{r) — acos'^ — -I- /3cos — , 
R R 
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with the same consequences concerning the explicit solvability and degeneracy. 

For the general case (|88| the same is true concerning the explicit solvability, 
however, as a rule, the degeneracy will be removed there even for the special 
case /3 = (i.e., /3 = 0, /3 = 0). 

Let us now go back to infinitesimal affinely-rigid bodies, thus, in addition to 
the gyroscopic degrees of freedom the deformative ones are taken into account. 
As mentioned, when the polar decomposition is used, then for the kinetic energy 
one obtains the formula ([50)1 with lo given by the expressions (pT|) and ([S^ . 
Just as in the gyroscopic case we concentrate on the phenomena in the two- 
dimensional spherical world. We use the same as previously parametrization of 
this world, i.e., (r, t^j))- or (i?, (/7)-coordinates. What concerns internal degrees of 
freedom we use coordinates -0, Cj where 



C/(^) = 



cos 
siu'f/; 



— sin?/^ 
cos V' 



and fj, A are factors of the polar decomposition. 
After some calculations we obtain that 



T = 




where the micromaterial coordinates are chosen (always possible) in such a way 
that the incrtial quadrupole is diagonal, i.e., 



J = 



Ji 






J2 



Obviously, lo is given by (|7^ . The first term of Tint above, i.e., the first term of 
([SD]) when n is specified to 2, is the centrifugal contribution; the coefficient at 
tt'^/2 is the dynamical deformation-dependent inertial momentum of the J7-rigid 
body (the material tensor AJ A in ([SO]) ). The second term is the Coriollis con- 
tribution to Tint, i.e., the coupling between angular and deformation velocities 
(Tr (AJ(dA/dt)cD) in ^^). And the third term of Tint is the kinetic energy of 
pure deformations <(Yx [j{dAI dtf) /2 in pi)) ). 

It is seen that even the purely kinetic term is rather complicated and there is 
no hope for integrability (or any kind of rigorous solvability) even if the potential 
term has a simple structure adapted to that of T. Obviously, the most natural 
potentials are those given by scalar invariants of the Green deformation tensor, 
i.e., analytically 



G 



e- 



-0 



7] 



C) 

,2 
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Obviously, in the case of some practical motivation some approximation or nu- 
merical procedures are possible, however, for the general J there is no hope for 
realistic potentials admitting integrability or some kind of analytical procedures. 

As usual, this is possible however for highly symmetric systems, when the 
internal inertia is isotropic, i.e.. 



J 



AB 



AB 



16 



AB 



This is explicitly seen when one uses the two-polar decomposition for internal 
degrees of freedom. 

For the general n we have the expression (1631) with given by the ex- 
pressions (|64|) . ((65|) . ((66|) . And there is also no hope for integrability or any 
kind of explicit analytical solvability. But again the special case n = 2, i.e., 
two-dimensional spherical world, is an exception. 

Parametrization of the S'^(0, i?)-world is again the same, i.e., (r, (p)- or ip)- 
variables. When expressed in terms of the two-polar decomposition, i.e.. 



ip = LDR-^ e GL(2,M), 
then p is parameterized by generalized coordinates a, (3, A, /i, where 



L(a) 



cos a 
sin a 



— sma 
cos a 



cos (3 — sin /? 
sin [3 cos (} 



DiX,fi) = 



A 
n 



In a partial analogy to ([72|) we obtain that 



^ r-idL 



-1 

1 



dt 



-1 

1 



where 



X = Xrl + Xdr 



da 



but d has no "drive" term, i.e., 



r dip 
R~dt 



da 



cosi? 



dip 



dt 



Remark: do not confuse this -d with the angle 

It turns out that to avoid some embarrassing cross-terms it is convenient to 
introduce the "mixed" coordinates 



1 



.:=-(A-M), 



1 



y 



V2 



a- f3. 



The inverse rules read that 
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The canonical momenta satisfy the contragradient rules 



^ {p\ -P^), py^ {p\ + Pf,) 



.y ^ 

Pi = l^{Pa+P(i) , P^"" 2 ~ ^'^^ 



and conversely, 



Pa^P-y+PS, Pa ^ Pi - PS- 

Remark: it is well known that the two-polar decomposition is not unique and 
has singular points [49]. Therefore, for n = 2 the (a, /3, A, /z)-parameterization is 
not a diffeomorphism between GL(2,IR) and x R^, i.e., the Cartesian product 
of the two-dimensional torus and the real plane. Rather the representation 
manifold x is plagued by a complicated system of identifications [JS] and 
only the resulting quotient space may be interpreted as a proper representation 
of our configuration space. What concerns the quantities 7, S the situation is 
even worse. They are "less proper" angular variables than a, (3 themselves are. 
The reason is that the integer-entries matrix transforming (a, (3) into (7, i5) has 
the determinant —2, thus, different from ±1. If we normalize it to the ±1 
determinant, then its entries will be no longer integer. Therefore, this matrix is 
not a well-defined automorphism either of the torus or the lattice 1? . The 
strange status of variables 7, 5 is not misleading if it is not forgotten. And 
of course, it is quite proper to say simply that to avoid some cross-terms in 
the kinetic energy expression we replace the angular velocities x,i? by their 
combinations, i.e.. 

Let us order our generalized coordinates q*, i = 1, 6, as follows: 
Then the kinetic energy will be written as follows: 



m dq' dqi 



where for the above ordering of variables the matrix [Gij] of the metric tensor 
G consists of three blocks subsequently placed along the diagonal (looking from 
the top to bottom): 

(i) the 1 X 1 block Mi, i.e.. 

Ml = [1] , 
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(a) the 3x3 block M2 given as follows: 



cos ^ 



—y cos ^ 

m ^ it 



T 9 
„ —X COS 
it m 



R m 



y cos 




{Hi) the 2x2 isotropic block M3, i.e., 



Mz = —h = 
m 



where, obviously, I2 denotes the 2x2 identity matrix. 

It is seen that now, when deformative degrees of freedom are taken into account, 
the previous special case / = niR^ does not lead to any essential simplification. 
One can easily show that 

det[Gij]=R' (^ly^Vsin^^, 

thus, the density of the Riemannian volume element equals 

,^ = i?(^)^.,|sin^. 

Explicitly, the block matrix [Gy] is given as follows: 

Ml 



Mo 



Ms 



The inverse contravariant tensor matrix [G*-'] is obviously given by 

" M-i 



M2-1 



M3-1 



where the inverse blocks have the forms: 



iii) M^' 



B? sin^ B? sin^ ^ 

__cosJ_ mill „i.„2T_ 

B? sin^ ^ / ^ ^ ;r^^^6 

__cosJ_ 1 2 r 

BP- sin^ i flS-Clg ^ 



cos 
■ fl2 sin" ^ 



o?Ctg o 



™ 1 4- 1 pfo-2 r 
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(lli) A/3-1 = fh = 

For potential systems with Lagrangians of the form 

L = T-V{q), 
the corresponding kinetic (geodetic) Hamiltonian equals 

2m 

and the full Hamiltonian is as follows: 

H = T + V{q). 

According to our convention of ordering coordinates q^, i = 1, 6, i.e., 

the corresponding conjugate momenta pi, i — 1,6, are denoted and ordered as 
follows: 

Pr,Pv,Pf,PS,Px,Py 

In certain expressions it is convenient to use the original momenta pa, p/s, p\, 
p^. First of all this concerns Pa,Pf3 because of their geometrical interpretation 
respectively as spin and vorticity [49] . 

For the above potential systems the reduced Hamilton- Jacobi equation has 
the form 

We do not write it down explicitly because of the complicated and rather obscure 
form of the resulting expression. The above block form is much more readable 
and lucid. 

Just as in the gyroscopic case we are dealing here with non-orthogonal coor- 
dinates (the 3x3 block M2) and it is not clear for us whether in some hypothetic 
orthonormal coordinates (they exist, of course) the system is separable. It is 
perhaps a little surprising that our kinetic Hamiltonian T has the separable 
structure. Of course, for the system with deformative degrees of freedom as 
above, the geodetic model is not physical because it admits unlimited expan- 
sion and contraction. Therefore, some potential must be assumed and this is 
just the problem, i.e., we could not determine a wide class of potentials com- 
patible with the separability in our non-orthogonal, but nevertheless natural, 
coordinates. Just as in the gyroscopic case we restrict ourselves to some special 
class of potentials, assuming in particular that all angles if, a, (3 (equivalently 
V3,7,(5) are cyclic variables. Therefore, the reduced action Sq is sought as a 
function linear in angular variables and separable, i.e., 

Soiq) = Sr{r) + S^ix) + Syiy) + iip + a,-f + CgS 

= SAr) + Sa,{x) + Sy{y)+£(p + Ca,a + C(3f3, (92) 



f 
f 
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where ^, C^, Cs, Con Cp are constants. The relationship between (7,(5) and 
{a, (3) imphes that 

i.e., ^ ^ 

^7=2 ^'^^ ' ~ 2 ~ ^'^'^ ' 

As Pa,P/3 corresponding respectively to the spin and vorticity in the flat-space 
theory [49], we shall also denote the constants Ca, Cp as s, j, thus, 

Soiq) = Sr{r) + Sa:{x) + Sy{y) + £ip + sa + j(3. 

Let us now quote the expressions for the action variables, i.e.. 





2ir 




Pada — 


Ca J da 


= 2'nCa = '2,7: s. 















Pfjdfi = 


Cp J dp 


= 2t:Cp = 27rj, 















p^dip = 


i j dip^ 












Let us assume that the potential energy separates explicitly with respect to 
acyclic variables, i.e., 

V{r,x,y) ^Vr{r) +V^(x) + Vy{y). 

Then the resulting Hamilton- Jacobi equation also separates and denoting the 
corresponding separation constants by Cr, Cx, Cy we have that 

C'r + Cx + Cy — E, 

and then 



+ Vr(r) = E-Cx-Cy, (93) 
+ Vxix) = Cx, (94) 

+ Vy{y) = Cy. (95) 



dSr^ 


2 




- JaCOS^)^ 


dr J 


)- 


87r2 


rni?2 sin^ 


hi 


^dSx 




(^ + Jpf 




y dx 




327r2/j;2 




( dSy 




( Jq — '//j)^ 




K dy 




327r2/?;2 
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Therefore, after expressing Ca — s, Cf3 — j, i with the help of J^, J/3, J^^ 

we obtain for the remaining action variables the expressions 



Due to the cyclic character of (p, a, the true dynamics is contained in (j93p . 
((94)) . ((95|) or on the level of action variables in jM]), Jill), ([Ml)- Obviously, the 
microdeformation dynamics is described by ([M]) . ([55)) or equivalently by ([W)) . 
(|55)) . Unlike this, (|M)) and (|55)) describe the dynamical influence of the "North 
Pole" (r — 0) and "South Pole" (r = -kK) on the translational motion. The 
dynamics of deformation invariants x, y is in many respects similar to the one 
in [37] (the combinations (Di - Di) /\/2 and (£»! + i:>2) /V2 in fS^). The r- 
geodetic model with K = is obviously well formulated. But in d'Alembert 
models the (x, j/)-geodetic case (T4 = 0, = 0) would be quite not physical 
because of admitting unlimited expansion and contraction of the body. This is 
not the case in affine models where the "elastic vibrations" may be encoded in 
the very kinetic energy. 

The procedure is now as follows. We assume some particular forms of 
the deformative potentials V^, Vy basing both on some physical arguments 
(taken, e.g., from elasticity theory) and on the possibility of explicit analyti- 
cal calculations (at least in terms of some familiar special functions) |37j. And 
later on, calculating ((96)) . (|97)) and inverting the resulting formulas, we ex- 
press Cx in terms of Jx-, J a, J (3 and Cy in terms of Jy, Ja, Jp- These expres- 
sions show some kind of "degeneracy" because Cx depends on Jq, Jp through 
the integer combination Ja + Jp, and similarly, Cy depends on them through 
Ja — Jp- However, this does not prejudice the true degeneracy of the Hamilto- 
nian. Substituting Cx [Jx, Ja + Jp) and Cy [Jy, Ja — Jp) to the main dynam- 
ical formula (|98)) and assuming again some functional form of Vr (e.g., just 

= in the translationally-free case), we in principle calculate Jr in terms of 
E, Ja, Jp, Jip, Jx, Jy And again, solving this expression with respect to E, one 
finds in principle that 



Without assuming the particular, both physical and analytically treatable, form 
of potentials one cannot a priori decide about problems concerning hyperinte- 
grability (degeneracy) of the above integrable models. 




(96) 



(97) 



(98) 





68 



Just as in the flat-space problems, there exist reasonably-looking models 
separable in other coordinates in the space of deformation invariants, moreover, 
separable simultaneously in several systems of coordinates in this space, thus 
probably degenerate (hyperintegrable) ones. 

An interesting class of separable models is obtained when one uses the polar 
coordinates {g, e) in the space of deformation invariants: 



X = Qsme, y=gcoss. 



(99) 



Then the kinetic energy 



TO 



'dt~dt 



with coordinates ordered like 

has the block matrix of the metric components 



[Gij] 



K3 



where 



Ki = [1], 



siv? ^ + —Q^ cos^ -5 —p^ s\T? £ cos ^ —p^ cos^ £ cos -5 



T 9 • 2 r 

— o sm ecos-k 

m'^ R 



— Q^ sin^e 



7" 9 9 r 

— Q cos e cos -k 

TO it 



— COS^ E 



-L 
^g^ 



(100) 



/ 
m 



1 
^2 



The inverse metric tensor underlying the kinetic part of the Hamiltonian is 
given as follows: 



where, obviously, 



K^' = [1], 



R^ sin^ 

cos 
~ _R2 sin^ 



I o2 



^Ctg2^ 



■ fl2 sin" ^ 



4ctg2^ 



^ + ^ctg2^ 



, (101) 
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Obviously, the purely kinetic (geodetic) Hamiltonian is given by 
where the canonical momenta are ordered as follows: 

iPl,P2,P3,P4,P5,P6) = {Pr,Pv,P'y,PS,Pg,Pe)- 

Strictly speaking, our "polar coordinates" q, e are rotated by 7r/2 with respect 
to the usual convention, but for some reasons it is convenient for us. 

Tensor densities built of the metric Gij are in these coordinates given by the 
expressions 

4 

4-2 2 -2 



G = det [Gij] = -R^ ( 3 ) sin^ e cos^ e sin^ , 



VG = R - 



/ \ 2 ^ 

— I Isinscosel sin — . 
mj R 

As previously, the system is separable (thus, completely integrable) for poten- 
tials independent of a, (3. Then the reduced action 5*0 is given by 

5o(g) = Sr{r) +1^ + C^a + Cpl3 + S^iQ) + 5e(e), 

where as previously £, Cq , C/3 are integration constants. Using again the terms 
of spin and vorticity s, j we can write that 

So{q) = Sr{r) +e^ + sa+ jP + 8^(9)+ Se{e). 

It is easily seen that such problems with cyclic variables ip, a, /3 are separable 
for deformation potentials of the form 

V,ig) + ^, (102) 

i.e., for the total potentials we have that 

V{r,Q,e) = Vr{r)+V,{Q) + ^. (103) 

Just as previously the action variables Ja, J/3, Jy, are given as follows: 

27r 



^ Pada = Oct J da = 2-iTCa = Stts, 


27r 

J0 = j> PffdP = Cp j dP = 27rC^ = 27rj, 


27r 

J^p = j> P(pdip = I J dip = 2n£. 
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Then Sg, S^, Sr satisfy the ordinary differential equations 



dSg 
dg 



+ 2Ig^iVgig)-C) 



-2IA. 



de 



Jl + 2 JaJ 13 cos2e + Jl 



1 



2m \ dr 



Att^ sin^ 2e 
dSr ^ ^ 



87r2mi?2 sin^ 



2IV,{e) - 21 A, 

E~C-Vr{r), 



where C and A are the separation constants. More precisely, 21 A is the constant 
value of the Hamilton- Jacobi term depending only on e, and {E — C) is the 
constant value of the r-tcrm. Therefore, the remaining action quantities are 
given by: 



J. = 



hl{A-V,{e))- 



Jl + 2 J^J (3 cos2e + Jj 



Att^ sin 2e 



de, 



/ 21 A 

l2l{C-v,{e}}-—de, 



\ 



2m{E - C ~Vr{r)) 



iyJip J a COS 

An^R^ sin^ 



-dr. 



(104) 
(105) 

(106) 



And again the calculation procedure is as follows. Some explicit forms of con- 
trolling potential functions Vg, Vr must be specified on the basis of mi- 
crophysical, phenomenological, or simply geometrical arguments. Later on one 
"calculates" the integral (|104p and expresses Jg as a function of A, Ja, Jp, i.e., 

Je = Je (A, Ja, Jp) ■ 

Solving this expression with respect to A, one obtains the formula for A as a 
function of three action variables: 

A = A{Ja,Jl3,Je)- 

Substituting this to pOSp and "calculating" the integral, one obtains Jg as a 
function of C, J a, Jp, Je, i-e., 

Jg = Jg (C, Ja, Jp, Je) ■ 

Solving this with respect to C, one can in principle express C as function of 
four action variables: 

C — C{Ja,Jp,Jg,Je)■ 
And finally, this expression is substituted to ()106p and then we obtain the 
formula for Jr as a function of Ja, Jp, Jg, Je, Jip, i-e., 



Jr — Jr{E, Ja, Jp, Jg, Je, Jip)- 
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Solving this with respect to E, one obtains the concluding formula expressing 
E through the (six) action variables: 

Unfortunately, unlike in the mechanics of material point moving in a central 
field, even the partial degeneracy cannot be concluded directly from the formulas 
(|104p . (|105p . (|106p without performing the integration process based on some 
fixed forms of the potential controlling functions V^, Vg, Vr- And just (and 
even more so) like in the mechanics of infinitesimal gyroscope, interaction with 
internal degrees of freedom removes the total degeneracy of the material-point 
Bertrand models ((55)) . (|7D)) . 

The potentials of the form (|102p . p03p are very convenient from the point 
of view of nonlinear macroscopic elasticity. Being compatible with the very 
nature of deformative degrees of freedom they are also interesting in the theory 
of infinitesimal objects. The following example is very instructive, moreover in 
macroscopic nonlinear elasticity in two dimensions it is almost canonical |49| : 

V=- - + -^,2^,^ _ + ^-^ , ,^>0. 107 

gi"'cos2e 2 2 / 

The first term prevents any kind of collapse of the two-dimensional body: to the 
point or to the straight line. The second term of the "harmonic oscillator" type 
prevents the unlimited expansion. The natural state A = /i = 1 (no deformation) 
minimizes the potential energy V] it is a stable equilibrium. Extension in one 
direction is accompanied by contraction in the orthogonal one. One can invent 
plenty of similar potentials just using the polar coordinates gi, e in the space of 
deformation invariants. The above one is particularly simple and suggestive; it 
is also well suited to the analytical procedure. 

5.2 Pseudospherical case 

Let us now consider the same problems on the pseudosphere, i.e., on the two- 
dimensional Lobatshevski space. We assume it to have the pseudoradius i?, 
therefore, the constant negative curvature TZ = —2/R^. The most natural real- 
ization of this space is the "upper" shell of the hyperboloid in R^, i.e., 

-x^ - y2 + = z > 0. 

We denote it as H^'^'^{0, R). Its metric tensor is obtained as the restriction 
(injection-pull-back) of the Minkowski metric: 

dx^ + dy^ — dz^ . 

Again we introduce the "polar" coordinates (r, ip) , where r is the geodetic dis- 
tance measured from the "North Pole" x = 0, y = 0, z = i?, and Lp is the 
"geographic longitude", i.e., the polar angle in the sense of the plane z = 0. 
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Remark: the "South Pole" x = 0, y = 0, z = —R is not interesting for 
us as it is placed on the "lower" shell of the hyperboloid, i.e., on the other 
connected component. H^'^'^ {0, R) is not compact and the radial variable r 
has the infinite range [0,oo]. We shall also use the pseudoangle = r/R. The 
parametric description of iJ^'^'+(0, R) C M.^ is given by 

x — Rshdcos(p, y = Rshi) simp, z = Rchd. 

Expressed in terms of coordinates (r, 1^9), the metric element of _ff^'^'+(0, R), in 
analogy to (|67|) . has the form 

R 

And in general, practically all pseudospherical formulas may be obtained from 
the spherical ones by substituting hyperbolic functions instead of the corre- 
sponding trigonometric ones. In certain expressions one must be however care- 
ful with signs. So, the translational kinetic energy, in analogy to (pB]) . has the 
form 

And again there are two Bertrand-type potentials [35l[40], i.e., 
{i) the "harmonic oscillator" -type potential: 

(ii) the "attractive Kepler-Coulomb" -type one: 

F(r) = --^cth^, a>0. (109) 
R R 

With these and only these potentials all bounded orbits are closed. And now the 
term "bounded" is essential because the "physical space" is now not compact. 
And indeed, there exist unbounded motions corresponding to energy values 
exceeding some thresholds. It is interesting that unlike in the spherical world, 
in Lobatshevski space the isotropic degenerate oscillator has an open subset of 
unbounded trajectories because the potential (|108p has a finite upper bound, 
i.e.. 

Sup V = 

For energy values above this threshold all trajectories are unbounded, the motion 
is infinite. Below this threshold all trajectories are not only bounded but also 
periodic. 

The existence of threshold in (|109[) is not surprising, it is like in the usual 
Kepler in M^. But the threshold for the isotropic degenerate oscillator is a very 
interesting feature of the Lobatshevski "world" . 
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Let us now consider an infinitesimal gyroscope. The success of ([7T|) suggests 
us to use the orthonormal basis 



or 



or in terms of components: 

Er = [1,0]^, 



1 



d 



1 



i?sh-^ dip i?sh 



" i?sh 



It is easy to see that this aholonomic basis is reahy orthonormal because 



9r 



9vr = 0, 



5r 



1, g^^ = i? sh 



R' 



We formally repeat all considerations of the spherical case, i.e., all formulas for 
Christoffel symbols are analogous, the trigonometric functions are replaced by 
the hyperbolic ones, etc. Preserving the same notation we obtain that 

dtp r dip d^j dip 
,^j+^,^. = _+eh-- = -+ch^-. 



T 



dip 


UJ 


'dt' 








^ Tint 


m 1 




'2 





dip 
R \~dt 



1 f dilj , r dip 
- — - + ch 

2 \dt R dt 
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I.e., 



T 



^tr + Tint 
TOi?2 f /d^ 

'dt 



V dt 



1 /dijj 

2 \~dt 



chi9 



dip 
'dt 



(111) 



When the generahzed coordinates {q^,q'^,q^^) are ordered as previously, i.e, 
{r,ip,ip), then the corresponding metric on the configuration space is given as 
follows: 

1 



i?2(sh' 







R. ^ mW^'^^ R) 



— ch- 



The kinetic energy is then given by the expression 







— ch— 

m R 



(112) 
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The contravariant metric [G"^] has the form 

10 " 

n '^^'g mil „iu2 r 

^ -WsiF^ T + IF'^^'^ "R . 

For potential systems the corresponding geodetic Hamihonian is written as fol- 
lows: 

2m 

where, obviously, the generalized momenta {pi,P2,P3) are ordered as follows: 
The weight-two density built of G has the form 

G' = det[Gy] sh^^, (113) 

and the weight-one density defining the volume element is given by 

VG = R\ —sh^. (114) 
V m R 

It is seen that the spherically very special case / = mR^ here, in the pseudo- 
spherical case also leads to some simplification of [Gij], but not so striking one 
as previously. This fact has deep geometric reasons which will be explained in 
the sequel. 

Namely, in the spherical space an essential point is the natural identifica- 
tion between the quotient manifold S0(3, M)/S0(2, R) and the spheres 5^(0, R), 
5*^(0, 1). And this has to do with the formal identification between two-dimen- 
sional rigid body moving over the spherical surface and the three-dimensional 
symmetrical top without translational degrees of freedom. The special case 
/ = mR^ corresponds to the spherical top. In general the kinetic energy is 
then invariant under S0(3,R) x S0(2,R). In the three-dimensional top analogy 
SO (3, R) is acting as left regular translations and SO (2, R) as right regular trans- 
lations corresponding to the group of rotations around the body-fixed z-axis. If 
/ = mi?^ we have the fuU invariance under S0(3,R) x S0(3,R). 

In the hyperbolic pseudospherical geometry the problem is isomorphic with 
the three-dimensional Lorentzian (Minkowskian) top on R^ . The rotation group 
S0(3,R) is replaced by the three-dimensional Lorentz group S0(l,2). And 
still an important role is played by S0(2,R) interpreted again as the group 
of usual rotations in Euclidean space of (x, ?/)-variables (thus, not affecting 
the z-direction) . The above kinetic energy (jllOp . (|llip is invariant under 
S0(l,2) X S0(2,R). But it is never invariant under S0(l,2) x S0(l,2), i.e., un- 
der left and right Lorentz regular translations in the S0(1, 2)-sense. The spher- 
ical special case / — mR^ does not help here. Indeed, the underlying metric G 
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(and the kinetic energy itself) is positively definite. But the doubly-invariant 
(S0(l,2) X S0(1, 2)-invariant) metric on S0(l,2), i.e., its Killing metric is not 
positively definite. Instead it has the normal-hyperbolic signature (-I--I — ). The 
reason is that it is semi-simple (even simple) non-compact group. This brings 
about the question about non-positive kinetic energies (metric tensors) on our 
configuration space. As the negative contribution to the Killing metric tensor 
on S0(l,2) comes from its compact subgroup S0(2,M) of (x, y)-rotations, i.e., 
from the gyroscopic degree of freedom in the language of 7J^'^'+(0, R), there is a 
natural suggestion to invert the sign of the gyroscopic contribution to ([75]) , ([7^ , 
i.e., to make it negative. One is naturally reluctant to indefinite kinetic energies 
but there are examples when they are just convenient and very useful as tools for 
describing some kinds of physical interactions [42l[44ll45ll46l[47ll48ll49l[50l[5T] . 
just encoding them even without any use of potentials. 

So, we can try to use, or at least mathematically analyze, the "Lorentz-type 
kinetic energies" Tl of the form 



mil 


' dr 








{{ 


2 





It 




2\dt R dt 

^f?+ch,?^ 
2 \ dt dt 
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Thus, it is so as if the extra rotation diminished effectively the kinetic energy 
of translational motion. If we write as usual that 

m dq' dqi 

TL-yLG.,(g) — — , 

then, with the same as previously convention concerning the ordering of coor- 
dinates (r, ip, -0), we have that 







i?2 (sJj2 r 







--ch-fe 







(compare this with (|112p ). And now, obviously, the remarkable simplification 
occurs in the very special case / = niR^ just as in the spherical symmetry. This 
has to do "as usual" with the enlarging of the symmetry group from S0(1, 2) x 
S0(2,M) to S0(l,2) X S0(l,2) (two additional parameters of symmetry). And 
namely, lG becomes then ^G, i.e.. 





-i?2ch^ 





-R^ch^ 
-R-" 



(compare this with ([82|) and notice the characteristic sign differences). 
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The weight- two scalar density buih of lG has the form 



det [iG„] =-i?2_sh2-, 
m R 



so it differs in sign from (jll3|) . This difference obviously does not influence the 
density of Riemannian volume element, i.e., 



det [iGi 



V m R 



exactly as in (|114p . This fact is of some interest. 

As usual, the geodetic Hamiltonian is given by the expression 

Tl = ^LG'^{q)p^Pj, 

where lG*-' are components of the contravariant inverse of LGij. They are given 
as follows: 

10 



[lG''] = 







1 



_cll^ 



ch-j 



rCth^ 



Obviously, if we use the above isomorphism between the two-dimensional top 
sliding over the Lobatshevski plane with the three-dimensional Lorentz top with- 
out translational motion in R^, then it is clear that ^G is, up to normalization, 
identical with the Killing metric tensor of SO (1,2). Let us quote some formu- 
las and concepts analogous to three-dimensional angular velocities, i.e., to (|76p . 
(|77|) . And then the kinetic energy will be expressed like in (|5T|) . 

First of all we parameterize S0(l,2) with the help of what we call the 
"pseudo-Euler angles" . So, let us write that 

S0(l,2) 3 = U,{ip)L,{d)U,{i>), 

where the meaning of Uz is like in (|75p and denotes some Lorentz transfor- 
mation in M^, namely, the "boost" along the x-axis, i.e.. 



1 

cM sM 
sM cW 



During the motion all these quantities are functions of time and we can calculate 
the corresponding Lie-algebraic element 
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i.e., the co-moving pseudo-angular velocity. After some calculations we obtain 
formulas analogous to ((771), ((78]), (HH), (HO]), and namely, 



A = Ai 



" 





" 




" 





1 " 
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-1 


" 








1 


+ A2 











+ A3 
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where 



Ai 



sh?? sin -0 



dip 



■ cos V' 



dt' 



A2 



A. 



cM 



shi? cos tfj- 
dip 



■ sin?/; 



dt '■ 



dt 

dip 

H 

(hp 

dt ' dt' 

Both the similarities and differences in comparison with the corresponding 
spherical formulas are easily seen. 

And now we can write two formulas analogous to 



T = — 



T = — 



^ A 



K 

K 
"2 



; A, 



^ A2 -:r A, 



(116) 
(117) 



where K > and / > 0. This is the symmetric S0(1, 2)-top in M^. The indef- 
inite expression (|117p is structurally suited to the normal-hyperbolic signature 
of S0(l,2). When K = I, then it becomes the spherical Lorentz top in M.^ in 
the indefinite version based on the Killing metric. 

It is easily seen that both expressions (|116|) and (|117p are invariant under 
S0(l,2) X S0(2,M), where S0(l,2) and S0(2,M) acts on S0(l,2) through re- 
spectively the left and right regular translations. The form (|117p with K = I 
is invariant under all regular translations (both left and right), i.e., under 
S0(l,2) X S0(l,2). And specifying K = mR^ in pTB]) and (fTTT]) . we obtain 
respectively piU)) and PT5)) . 

Let us now turn to the Hamilton-Jacobi equation and action-angle variables. 
Just as previously, ip, ?/> are cyclic variables in the kinetic energy term and we 
assume that the same is true for the total Hamiltonian, i.e., that the potential 
energy depends only on the variable r. Because of the non-diagonal structure of 
[Gy ] and [lG^] in the natural variables (r, (p, ?/;), it is rather difficult to decide 
a priori what would be (if any!) the form of the general separable potential. So, 
just as in ((85|) the reduced action is sought (with the same meaning of all 
symbols) in the form 

5*0 (r, ip, -0; E, £, s) — Sr{r, E) + £ip + sip. 

Just as in ([86|) the genuine dynamics reduces to the radial function which sat- 
isfies the equation 

2 /n , 



dSr 

dr 



2m {E-V{r)) - 



i?2sh2^ 



± 



TO 



(118) 
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The ± signs in (|118p refer respectively to the models (jllOp . (jllSp . Just as in 
the spherical case Pi^, are constants of motion and their constant values on 
fixed trajectories equal 

dSo . dSo 
And, of course, their action variables are given as follows: 

27r 

Jip = <f Pipdip — / id(p — 2tt£, 





In 



p^dip =^ J sdi/j ~ 2tts. 



Substituting this to the radial action expression and to ()118p we obtain for 

Jr — <j) Prdr 

the following formula: 



— J^ch-^) m J, 



2 



2m{E - V{r)) - ^ IZI^ii^ ± '^i^dr (119) 



/ 47r 



(with the above-mentioned meaning of the ± signs). 

From now on the procedure is just as previously. After substituting here 
some particular form of V we perform the above contour integration and (at 
least in principle) determine the dependence of Jr on E, J^, J^, i.e., 

Jr — Jri^J^^ J ip ^ J^^' 

Solving this equation with respect to E', one obtains the final expression for 
Hamiltonian: 

E = H{Jr, Jip, Jifj)- 

Just as in the spherical case no degeneracy structure may be directly deduced 
from (|119[) even in the special case of the minus sign and / — mR^ . 

An important point is that, because of the uncompactness of the Lobat- 
shevski space, the geodetic motion {V — 0) is unbounded, and then obviously 
the action-angle formalism becomes meaningless. Nevertheless, even for un- 
bounded motion the isotropic models with V depending only on r may be effec- 
tively studied with the use of explicitly known constants of motion and reduced 
to first-order ordinary differential equations. 

Let us now consider a deformable top moving in Lobatshevski space. All 
symbols concerning internal degrees of freedom arc just those used in spherical 
geometry. The metric tensor G underlying the kinetic energy expression, i.e., 

m dq' dqi 
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has the form analogous to the spherical case with the trigonometric functions 
simply replaced by the hyperbolic ones without any change of sign. Thus, the 
matrix [Gij] consists of three blocks Mi,M2,Ms, where 



Ml = [1], 



R m \ ^ J R m 



^x^ch^ 

m it 



-2/^ch^ 

m. «^ I 



M3 = —h 

m 







m. ^ 



When deformations are admitted, then there is no particular geometric mo- 
tivation for discussing the "Lorentzian" model with the indefinite kinetic energy. 
This is nevertheless possible and formally simply consists in replacing / by — 7 
in all formulas; in other words nc\a;ative "incrtial moments" arc admitted. 

For the inverse contravariant metric [G*^] underlying the geodetic Hamilto- 
nian, i.e., 

we have the block structure also quite analogous to the spherical formulas: 

Mfi = [1], 



ch 



_ch^ 



ch-^ 



ch 



ft 



^ ~1~ cth 



I y2 



+ ;^Cth^ 



Mo- 



rn . 



f 



f 



.3 -J^2- 

The corresponding scalar G-densities are given as follows: 



G = det [Gij] 
x/\G\ 



R 



m 
I 
m 



\xy\ sh 



R' 



Again we assume that the potential energy docs not depend on the angles 
If, a, P, i.e., they are cyclic variables for the total Hamiltonian. So, just as 
previously, the complete integral of 



E 
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will be sought as a linear function of all angles, i.e., 

Soiq) = Sr{r) + S^ix) + Syiy) + iip + a,-f + CsS 
= Sr{r) + S^{x) + Sy{y)+£ip + Caa + Cf}f3, 

with the same as previously relationships between constants. Coming back to 
the more familiar spin and vorticity symbols s,j, we have that C'a — s and 
= j. The action variables corresponding to Jq,, J^, are respectively given 
by the same formulas like in the spherical geometry: 



Jb = 





2-K 


Pada — 


Ca J da 











Ppdf3 = 


j dp . 











p^dLp — 


i J dip = 








Similarly, when the {x, y)-deformation invariants are used, the most natural 
separable potentials have the explicitly separated form: 

Vir,x,y) ^Vrir) +V,ix) + Vy{y). 

Then, just as in the spherical case, denoting the separation constants by Cr, 
Cx, and Cy, we have that 

Cr + Cx + Cy — E 

and the expressions for Jr,Jx, Jy just analogous to ([M)) . ([W]) . (1551) . namely, 



Jx ^ (P Pxdx = \l2I{Cx - Vx{x)) - ^"^l^J^}^ dx. 



Jy = ^p^dy=<^^\2I{Cy~Vy{y))^^-^^^dv, 



r \2 



J, ^ <j> prdr = ^ y 2m{E - Cx - Cy - Vr{r)) - ' J^2j^2°fa2i 

The procedure of eliminating constants Cx, Cy and expressing E through the 
action variables, i.e., 

E — l~L{^Jr : Jip , J a, <-^/3; Jx: '-^y) ^ 

proceeds exactly as in the spherical case. 
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Exactly as in the theory of deformable gyroscope in the spherical space 
it is convenient and practically useful to parameterize deformation invariants 
with the use of polar variables g, e (see ([M)) '). The only formal difference is 
that the trigonometric functions oi r/R (but not those of e!) are replaced by 
the hyperbolic ones without the change of sign, thus, (jlOOp and (|10ip take on 
respectively the forms 



R^sh' ^ + -L p^ch^ ^ -Lg^ sin^ ech^ ^ g' cos' ech^ 



K- 



— sin^ £ch-^ 

i^e' cos^ech^ 
1 



/ 2 ■ 2 

— g sm e 



_chj_ 



irth^^ 



2 2 

cos^e 

__chj_ 



rCth^ 



rCth^ 



/ cos^ e 



-^cth^ 



And obviously the other blocks are identical with those from the spherical ge- 
ometry. The metrical scalar densities are given as follows: 



4-2 2 1,1 

g sm ecos e sh 



i?' 



G = det [G,j ] = i?^ 

^|g[ = ^ ( £: ) Isinecosej sh— . 



There are also no essential changes with the integration procedure, reasonable 
deformative potentials models, and the action-angle variables. Just the trigono- 
metric functions oi r/R replaced by the hyperbolic ones. 



5.3 Toroidal case 

Let us now consider another interesting two-dimensional model, namely, the 
motion of structured material points over the toroidal manifold, so this time 
over the Riemann space of not constant curvature. Obviously, we do not mean 
the torus M^/Z^, i.e., the quotient of with respect to the "crystal lattice" 
1?; this would be flat, although topologically not trivial, space. Our torus is a 
"tyre" , or rather "inner tube" , injected into and endowed with the metric 
tensor induced from 'E? (from its flat Euclidean structure) . 

Let R denote the "small" radius, i.e., the radius of circles obtained from 
transversal orthogonal cross-sections of the tube by planes. The "large" radius, 
i.e., the radius of the centrally placed inner-tube circle will be denoted by L. 
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The parametric description of the torus is given as follows: 



X = {L + R cos 1?) cos ip, 
y = {L + R cos 1?) sin tp, 
z = iJsint?, 

where (p is the angle measured along the central inner circle ( "along- tyre" angle) 
and 1? denotes the polar angle of cross-section circles ( "around- tyre" angle) 
measured from the external circle (maximally remote from the "tyre" centre). 
Both ip and i3 run over the usual angular range [0, 2tt] with its obvious non- 
uniqueness. One can easily show that such a torus, denoted by T^{0,L,R), 
consists of points {x, y, z) € satisfying the fourth-degree algebraic equation 

(x^ +y^ + z^ + - R^'f - AL^{x^ + y^) = Q. 

Therefore, it is an algebraic manifold of fourth-degree because it is impossible 
to lower the degree of the polynomial on the left-hand side. 

After the easy calculation one obtains that the metric element induced on 
the surface r^(0, L, R) has the form 

ds'^ = R^M'^ + {L + Rcos'df dip'^. 

Introducing as previously the geodetic length along the " around- tyre" circles, 
i.e., 

r = Rd& [0, 2ttR\ , 

we can write that 



ds^ = dr^ + {l + R cos ^ dip^ . 



Therefore, the kinetic energy of the structure-less material point is given by the 
expression 





) +{L + Rcos^f 



In analogy to spherical and pseudo-spherical formulas it is also convenient 
to use the form 



^ + - -H cost? 





It is seen that again tp is the cyclic variable in T^,.. The along-(/9 rotations, i.e., 
ip ^ ip + a, are isometrics, thus, d/dip is the Killing vector field. 
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As previously, the coordinates on T^{0,L,R) will be ordered as (i?, 97) or 
{r,ip), thus, 



[9ij] 



i?2 
(i + i?C0St?)2 



1 

W 




{L+R COS 

and the density of Riemannian volume (area) is given as follows: 

y^\ = R{L + R cos ^) . 
Separable potentials have the form 



F(i?,(p) = K?(i?) + 



(L + i? cos 



The simplest and geometrically interesting is the class of models invariant under 
(^-rotations, i.e., = 0. The angle (p is then a cyclic variable. Unfortunately, 
even the simplest case, i.e., the geodetic model V = 0, is technically rather 
complicated and leads to elliptic integrals (after an appropriate change of coor- 
dinates). 

One can easily show that the only not vanishing Christoffel symbols are as 
follows: 



-pip 



Rsin'd 
"L + iJcosi?' 



— + cos 1? I sin a?. 
R 



As expected, the most natural choice of the auxiliary field of frames E is one 
suited to coordinate lines, i.e.. 



^""-Rd^-R^"' 



E,r. 



1 



d 



1 



L + Rcos'd dtp L + Rcos'd 



It is evidently orthogonal because 



g^dip = 9<P'0 = 0, g^^ = R^, g^^ = {L + Rcos'&Y . 

Let us now take into account internal degrees of freedom. We begin with the 
infinitesimal gyroscope. The angular velocities are then analytically expressed 
by the expressions 



-^ + sm'd-f. 
dt dt 



Therefore, the kinetic energy of the infinitesimal rotator is given as follows: 

T = Ttr + Tint 



m 



d-d 



dt 



- +(L + i?cos^?)^ 



d(p 



dt 



I (d^ . d^' 
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If we write it in the previous form, i.e., 



ltt~dt' 



with the ordering 
then 



{q',q\q') = {^,^,i^), 



E? 



{L + Rcos'df + j-sm^ d ^sini? 







-^-sini? 



The inverse metric underlying the geodetic Hamiltonian is given by the expres- 
sion 

" ^ 



1 

BP- 



n 1 

" (L+JJcosi?)^ 







sin xf 
' (L+Jicosi?)^ 



(L+Rcos^y I ^ (L+Rcos-dy 

Up to the /-dependent normahzation, the Riemannian density function is 
like that for the structure-less material point, i.e., 



IGI = \ — RiL + Rcos'd) . 
V m 

Just as in the spherical and pseudospherical spaces, the metric G is not diag- 
onal in natural coordinates, therefore, without rather complicated analysis it 
would be difficult to fix the class of potentials compatible with the separation of 
variables procedure. And again evidently integrable models are isotropic ones 
when the variables ip, are cyclic, i.e., 

H = T + V{^) = ^G'^PiPj + 
Obviously, the phase space coordinates are ordered in the usual way, i.e.. 

For the action variables J^, we obtain that 

271- 





^ p^dip = j . 



sdtp = 2tts, 



(120) 
(121) 
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with the same meaning of symbols as previously and 

The true dynamics is contained in J§ which after the substitution of (jl20p and 
(fnT|) into 

becomes as follows: 



^ ^' A^-^ (L + Rcosdf I 47r2 

Assuming some particular shape of V{'d)^ "performing" the integration, and 
solving the result with respect to E, one obtains the final formula: 

E — Ti (Jij, J^p, J^) . 

No degeneracy is a priori seen without performing the integration, even for the 
geodetic models, i.e., when V = Q. An elementary integrability (if any!) may 
be suspected, of course, only for ¥{•&) built in a simple way from sini^, cosi9- 
expressions. 

Let us now quote some formulas concerning an infinitesimal deformable top 
moving over the toroidal surface. We restrict ourselves to the doubly isotropic 
dynamical models of internal degrees of freedom and use the two-polar decom- 
position of the matrix Lp. All symbols are just as in the spherical and pseudo- 
spherical case, in particular, generalized coordinates g' , i = 1, 6, are ordered as 
follows: 

Angular velocities x, 9 of the two-polar decomposition are given by the expres- 
sions 

da . diy9 dj3 
X = X.i + Xd. = ^+smz?-, 9=-, 

where obviously a, (3 are the primary angles of the decomposition, i.e., 
a = i (7 + 5) , /3 = i (7 - J) . 
The general formula implies then that 

r = r,, + T,„, = -G.,(g)— — , 

where with the above ordering of generalized coordinates the matrix [Gij] con- 
sists of three distinguished blocks Mi, M2, M3, and namely, 



86 



Ah 







—x^ sin?? 

m 

— sin^ -d 

m ^ 







X„2 



m 



The Riemann density has the form 

yjG|=i?(^i-^ |a;2;|(L + i?cos?9). 
The inverse matrix [G*-'] consists of the blocks 



i?2 

sin t9 



11? 



(L+Rcos^)'-' 

sin 1? 

JL+Rcosly^ 

sin 1? 

(L+flcosiS)2 



mil sin^ 1? 

I ^ jL+RcosWj^ 

sin^ i5 



(L+ficosi?)^ 



sin^ 1? 
(L+i?cosi9)^ 



(L+_RcosiJ)2 

f 



m 1 
/ y2 



sin^ i9 



(L+flcosi9)2 







As in all previous examples, the G-non-orthogonality of the natural coordinates 
is a serious difficulty in analysis of integrability problems. And just as previously, 
one can do something when the angles ip, a, f3 are cyclic variables, i.e., when 
V depends only on "i? and deformations invariants x, y. Namely, the system is 
separable for Hamiltonians of the form 



H ^T+V4^) + V^ix)+Vy{y). 



(122) 



When Vx , Vy are chosen so that to prevent an unlimited extension or contraction 
of the body, then they are capable to encode something like the dynamics of 
nonlinear elastic vibrations. This is even much more true when we use the polar 
coordinates {g,s), i.e., (I99p . on the plane of deformation invariants: 



X = gsins, 



y = gcoss. 



In these coordinates the matrix [Gij] consists again of three blocks Ki, K2, K3, 
where Ki = Mi, K2 is just M2 with formally substituted ([M]) . and if 3 has the 
form 



± 







m 



/ 

TO 



1 

o2 
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Just as in the spherical and pseudospherical case the problem is separable for 
deformation potentials (|102p : 

Q 

i.e., for the total potentials 

In particular, (|107p is very interesting from the point of view of nonlinear elas- 
ticity. 

The general procedure of constructing the action J-variables and discussing 
degeneracy problems is just like in the spherical and pseudospherical cases. Let 
us write down some explicit formulas. 

According to our assumption about independence of V of the angles Lp, a, (3 
we have the following forms of the reduced action 5*0: 

5*0 ^ S^['d) + S^{x) + Sy{y)+tif + sa + jP 
= Sii{d) + Se{Q) + S,{e)+l^ + sa + j(}, 

depending on whether we use the (x, y)- or (gi, e)-variables in the plane of de- 
formations invariants. Therefore, just as previously. 





27r 




Pada = 


J sda = 


27rs, 















P lid 13 = 


J jdP = 


27rj, 











27r 




p^difi = 


j idip = 


27r£, 



where s, j, £ denote respectively the constant fixed values of spin, vorticity, and 
"orbital" angular momenta. 

If Hamiltonian has the form (|122p , then the separation of variables procedure 
tells us that 

J. ^ / W 2/ (C. - - ^''l^J^} dx, (123) 



Jy^f pliCy - Vy{y)) i:^^-^ rfj/ , (124) 

Ji)^R<b j2m{E-C,-Cy~ Kj(i9)) - ('^'p " sm^)' 
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And again the procedure is like previously. When some explicit forms oiVx, Vy 
are assumed, the integrals (|123p . (|124p are "calculated", and then Cx, C, 
are expressed as functions of Jx, Ja, Jp, i-e.. 



Substituting this to (|125p and "calculating" the integral, we express J§ through 
E, Ja, J 13, Jtp, Jx, Jy And finally, solving this expression with respect to E, 
we determine E as a function of action variables, i.e.. 



And unfortunately again without explicit calculations nothing may be decided 
a priori concerning even some partial degeneracy. 

6 Three-dimensional problems 

Now we are going to discuss the motion of structured material points in three 
dimensions. We concentrate on the very interesting special case of the Einstein 
Universe, i.e., compact constant-curvature space metrically diffeomorphic with 
the three-dimensional sphere S^{0,R) C K*. This three-dimensional model is 
quite different from the previously studied two-dimensional ones. The point is 
that the spheres S^{0,R) are parallehzable, moreover, 5*^(0, 1) may be identi- 
fied in a natural way with the group SU(2). This means that the most natural 
choice of the auxiliary reference frame E is that of basic generators of left or 
right regular translations on SU(2), i.e., respectively, the basic system of right- 
or left-invariant vector fields on SU(2). This frame may be easily chosen or- 
thonormal. And all this means that the principal bundle of orthonormal frames 
over SU(2), FM{S\J{2), g) may be naturally identified with the Cartesian prod- 
uct SU(2) X S0(3,M). The metric g here is obtained by the natural restriction 
of the Euclidean metric in R'* to 5*^(0, R) and in the SU(2)-terms it is identical 
up to a constant factor with the Killing-Cartan metric on SU(2). The latter is 
invariant under left and right regular translations on the group SU(2). 

Before going any further we return temporarily to the more general n- 
dimensional formulation. This makes our description more lucid. So, we return 
to the general formula ([5T|l for the co-moving affine velocity VI, i.e, 






" A = F-l- DCf eV + If 



dt 



(126) 



Let us remind that 17 is defined by the formula 



DcA 
Dt 



I.e. 
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When the internal motion is gyroscopic, i.e., E SO(n,M), then obviously is 
skew-symmetric in the Kronecker-delta sense, i.e.. 

In the above formula (|126|) denote the co-moving components of the trans- 
lational velocity, i.e., 

dt 

And we shall use the abbreviations ([FT]) , i.e., 

f^d/fi = v-^''fT^dcv''b^''eV'', (127) 

i'rl B — </5 C T — , "rl B — T)— B, U^oJ 

dt dt 



thus, 



Q = r^dr + fill = ^di + "yS ^^rl(p. 



If, as we always assume, the auxiliary frame E is g-orthonormal and T is the 
(7-Levi-Civita connection, or more generally it is some Riemann-Cartan connec- 
tion, then the aholonomic components T^bc stre Kronecker-skew-symmetric in 
the first pair of indices, i.e., 

T^A A r cALt^K 

i BC = B C = —ObkO i LC- 

It is not true for the holonomic components T^jk because they are not g-skew- 
symmetric, 

bk- 



Let us now assume that the "legs" Ea of the auxiliary field of frames E form 
a Lie algebra, i.e., 

[Ea,Eb]--C''abEk, 

where ab are some constants, just the structure constants of the algebra 
with respect to the basis (. . . , Ea, ■ ■ ■)■ Lie bracket of vector fields in the above 
formula is meant in the convention 

i.e., as the commutator of X, Y when the vector fields are identified with first- 
order differential operators: 

xf - X7^,. 

If some point xq € M is fixed, then the manifold M may be identified with 
some Lie group and Xq with its identity. Without distinguished xq the manifold 
M is a homogeneous space of the corresponding Lie group with trivial isotropy 
groups. In the Abelian case with the M"-topology of AI this is just the difference 
between the linear and affine spaces. 
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The algebraic Killing metric of the Lie algebra with structure constants 
LM has with respect to the basis (. . . , Ea, . . .) the components 

lAB — laC^ KB- 

The corresponding Killing tensor field on M is given in local coordinates as 
follows: 

l^j =-fABE'^^E^j, 

i.e., in the absolute tensor notation we have that 

7 - IabE^ ® E^. 

Obviously, it is not singular if and only if the Lie algebra given by structure 
constants C is semi-simple. It has the elliptic signature, i.e., is essentially Rie- 
mannian when the underlying Lie group is compact. Then strictly speaking the 
metric is negatively definite, but changing its over-all sign we obtain the usual 
positive metric. The inverse algebraic metric satisfies the condition 

1 ICB 

thus, the contravariant metric field on M is given as follows: 

= E'aE^bI^'', 

or in the index-free form: 

^ = -i^^EA®EB. 

While ^ijE^ aE^ B = Iab — const, then the metric tensors 7, g on M define 
practically "the same" geometries. And from now on we assume that 7 and g 
just coincide at least up to a constant multiplier. Let us mention by the way 
that this multiplier does not affect the Levi-Civita connection. 

Let r[£'] denote the teleparallelism affine connection (|45)) . (|46)) built of E 
and S[E] is its torsion tensor (gT]). It follows from (gH), that S[E]^bc, i-e., 
the aholonomic coefficients of S[E] with respect to E, are given as follows: 

S [E]^ BC = -^C^BC- 

We have just assumed that the metric tensors 7, g differ by an overall constant 
multiplier A, thus, 

jab = A^As, li] = A^y. (129) 

Their common Levi-Civita connection will be denoted by {} or componentwisely 

by Tc !'■ -'^'^'^'^^'diirg to the assumed relationships between g, E, 7, the metric 

tensors 7, g are parallel under the connection T[E] (i.e., the teleparallelism 
connection is metrical with respect to 7, g). Therefore, following p3p . 

r [Ef jk = s > + s'jk + + Skj\ 
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i.e., in the aholonomic terms: 



r [E]'^ LM - I j^j^j I + -jC^ LM + -jCLM^ - -jCu^ L- 

The shift of "small" holonomic indices is meant here in the sense of gij , whereas 
the shift of aholonomic "capital" indices is meant in the Kronecker-delta sense. 
But the i?-aholonomic components of the £'-teleparallelism connection evidently 
vanish, i.e., r[EY^ lm — 0, so we obtain that 

E I 1/-. A' , 1/-. K 

LM — T^^LM + -^^M L- 



LM j 2 2 2 

In semi-simple Lie algebras the structure constants are skew-symmetric in the 
first pair of indices with respect to the Killing metric. As we have assumed the 
latter to be proportional to the Kronecker delta, then we finally obtain that 

'^M L — —7:^ LM- 



LM J 2 2 
So, finally ([T?7)) . ([T^ imply that 

= ^ eV +V C ■ 

It turns out that in certain formulas for kinetic energy it may be convenient to 
use other representations for the relative (internal) aholonomic velocity. Namely, 
from the point of view of GL(n, R)-objects firi is the "co-moving" affine/angular 
velocity. One can also use the "spatial" representation in R", i.e., 

^ d(p^c ic ,A o c ,„-l-D 
"rl B = — T— 'Z' B = ^ C"rl D'fi B- 

at 

There are also some additional possibilities concerning the M-representation, 
i.e., literally understood spatial representation. Obviously, the most geometric 
one is that based on the e- injection, i.e.. 



o* — 0^4 „B _ 



Dt 

o i o ^ ^B ,i A ,A ^B 

ilri J = e A"ri se J, ip J = e BS j, 



and so on. But one can also use the i?-injection of M" into the tangent spaces 
of M, i.e., 



r!ri', = i;^Af^ri i, v' = eW', ^', = e'av''bE",, 

and so on. 
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Let us now go back to the three-dimensional simple compact groups SU(2), 
S0(3,K). In the standard representation structure constants are given by the 
three-dimensional Ricci symbol: 



LM 



LM, 



Cklm = £klm 



(again the shift of indices is meant in the Kronecker sense). 

Some analytical expressions must be quoted now. Lie algebra of SU(2) con- 
sists of trace-less anti-hermitian 2x2 matrices and it is customary to use the 
standard basis a G SU(2)', j — 1, 2, 3, where 



1 

Zl 



and fjB denote the Pauli matrices: 





■ 


1 ■ 




■ -i' 




■ 1 " 




1 





0-2 = 


i 




-1 



Using canonical coordinates of the first kind k E 
for the SU(2)-elements 



, (A:) = exp (A; as) = /2 cos — — icrs— sin —, 



, we have the expressions 
(130) 



where I2 



1 
1 



is the 2x2 identity matrix and 



is the Euclidean length of k. The direction k/k G 5*^(0, 1) of k runs over the 
total indicated range, whereas the modulus k changes between and 2tt. All 
vectors k with k = 2tt refer to the same group element —I2 G SU(2). For the 
values of k larger than 27r the formula (|130p would repeat the former elements 
of SU(2). In this way A: = and k = 211 are singularities of the above coordinate 
system. Obviously, aj obey the standard commutation rules, i.e.. 



.M 



The Lie algebra of SO (3, R), i.e., SO (3, R)', consists of real 3x3 skew-symmetric 
matrices. The standard basis is given by Ak, K = 1,2, 3, where 



^ L 
'£k M, 



and the trivial (purely cosmetic) Kronecker shift of indices is meant. With this 
convention wc have the same commutation rules, i.e.. 



[Aj,AK]=e''jKA 



M- 
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The exponential mapping 

R (fc) = exp {k^Aj) (131) 
leads to the explicit form of finite rotation matrices R {li) : 

R(k) -x — cos k -x + {\ — cos /c) ■x \ ^ + sin/c ^ x 5, 
^ ' \k J k k 

with the standard notation a ■ b and a x 6 respectively for the scalar and vector 
product in M?. 

Now k runs over the range [0, tt] and antipodal points on the coordinate 
sphere A; = tt in are pairwise identified because 

R{Trn) — R{—Trn) 

for any versor n g M'^, n ■ n = 1. In this way k is the usual rotation vector, 
k is the rotation angle (in radians), and fl = fc/fc is the rotation axis in the 
right-hand-screw sense, R(k)k = k. In particular, we have the intuitive formula 

R(k)x ^x + kxx+ ^k X (k xx) + ■ ■ ■ + -^k x (k x ■■■ (k x x) ...) + ■■ ■ , 
2! ?i! 

where the term with the l/7i!-factor contains n copies of k. The natural homo- 
morphism 

pr : SU(2) ^ S0(3,M) 

is given hy v t—^ R, where 

vu{k)v~'^ = u[Rk). 

Therefore, 



pr 



-^{R{k)) = i^u{k),u[~{2^ ^ k)^"^ = {u(k),-u{k)} . 



The formula (|130p injects the group SU(2) into the real four-dimensional 
space spanned by matrices /2, —icrs, B — 1,2,3. Obviously, this space may 
be identified with R^. It is seen that the sum of squared coefficients at the 
mentioned matrices equals one. This is just the mentioned identification between 
SU(2) and three-dimensional unit sphere 5''^(0, 1) C K"*. 

The Lie- algebraic Killing metric built of the structure constants bc has 
the components 

lAB = -2(5ab, 

therefore, in (|129p we have that 

A = -2. 

The metric field on SU(2) or S0(3,M) is analytically given as follows: 

5,^__sm -J,,.+ M-_sm--j--, (132) 
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i.e., the corresponding arc element has the form 

ds^ = dfc^ + 4 sin^ ^ [d'd'^ + sin^ Mip^) , (133) 

where k, (p are spherical variables in the space of rotation vectors k. The 
angular variables i?, ip parameterize the manifold of rotation versors n{d, (p) = 
k/k and we can write that 

k 

ds^ = dk^ + 4 sin^ —dn ■ dn, 

or using more sophisticated tensorial terms: 

k 

g = dkiSidk + 4: sin^ -SAsdn^ dn^ . 

The reference frame E may be chosen either as ''E or ^E respectively con- 
sisting of vector fields generating left or right translations on SU(2), S0(3,M). 
The vector fields ^Ea, ^Ea arc respectively right- and left-invariant. Let ^E, "^E 
be their dual co-frames. The covector fields (differential one-forms) ''E^, ^E^ 
are respectively right- and left-invariant. 

After some calculations one obtains the analytical expressions in terms of 
coordinates k: 

k ^ k d k ^ k\kAk-^ d 1 j,m 



= r'^-20k^ + v-r'^-2)TTdP + r--'8P^ 

k k d / k k\kAk-' d 1 j,m 9 
'Ea = i;ctg--^ + {l--ctg-]— — — J- -SAM'k^ 



2 ^2dk^ \ 2 ^2J k k 2 ^™ 8^' 

sin fc „ 4 sink\ k^ kB B ^ ■ 2 f^- a , b „n 
^E^ = ——dk^+il — ] —^dk^ + —rum'^ - s^Bck dk^, 



k V k J k k fc2 2 

,^A ^ ^^dk^+(i_^j^\^^dk'' -^sin^-e^Bck'^dk^ 
k \ k I k k k'^ 2 



The vector fields 



DA~^EA-''EA=eAB^k'' ^ 



are infinitesimal generators of inner automorphisms, i.e., 

U VUV~^. 

It is convenient to separate explicitly the "radial" variable k and the angles 
1?, (f. To do this in a symmetric way, we use differential operators (vector fields) 
Da which depend only on t?, (p and act only on these variables, i.e., 

DAf = 0, 

if / is a function of the variable k. Similarly, the redundant system of quantities 
n = k/k are used. Obviously, n"^ depend only on i?, (p, but not on k. 
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Then one can show that 



d 1 k 1 
^Ea = UAgj^ - -ctg- EABcn^D^ + -Da, 

d 1 k 1 
''Ea = riA^ - 2*^*82 ^"^Bcn^D^ - -^Da, 

k 

''E^ = n^dfc + 2sin^ - e^-^'^nsdnc + sinfc dn^, 
k 

''E^ = n^dfc - 2sin^ - e^-^'^nsdnc + sinA; dn^, 



or in a brief, Httle symbohc way 



_d 1 fc _ — 1 — 

_ 9 1 fc _ _ 1_ 

E — n— eta— n X Z? £), 

(9fc 2 ^2 2 ' 

~ fc 
'i? = Tidfc + 2 sin^ — n X c?rT + sin fc rfn, 

~ fc 
''i? = ridk — 2 sin^ —nX(]m,~\' sin fc oBi. 



Obviously, 



{dk,DA)^DAk = 0, (dk,-^^\, 

d \ ^72 

afc / ^ "SF " ^' ^'^^ ^ ^B"^^ = eabc"-'^- 

Commutation relations between the above vector fields read as follows: 
[^Ea,''Eb] — ~£ab'" ^Ec, YEa,^Eb\ — Sab" ^Ec-, 
[EaJEb] -0. 

The vanishing of the last Lie bracket is due to the fact that the left and right 
group translations mutually commute. 

The metric tensor g (see (|132p . (IISSI) ). i.e., the (— 2)-multiplc of the Killing 
tensor, may be written down as follows: 

g = SabE^ ® 'E^ = 6ab'^E^ ® 

Similarly, the contravariant inverse metric 

fc2 / fc2 \ 



4 sin^ I V ^ ^^^^ t 



may be expressed as follows: 

g = S^^'Ea (E> 'Eb = S^^'-Ea ® "^Eb, 
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or using the k, n- variables: 

~ d d 1 ^ATj ^ ^ 

5 = €3 ^ + —^S^^Da ® Db. 
ok dk 4 sm^ f 

Let us now pass over from the standard unit sphere S'^(0, 1) ~ SU(2) to the 
sphere S^{Q,R) C M"^ of the radius R. We introduce the new "radial" variable 
r = Rk/2, i.e., k = 2r/R. Obviously, the range of r is [0,7ri?]. The coordinate 
vector k is replaced by f with the length r and such that 

r k _ 
r k 

Expressions for the auxiliary reference frames E and co-frames E take on the 
forms 

R R 

Ea = E{R)a, ^Ea = —'^E{R)a, 

= ^^E{R)^, "-E"^ = ^'■E{R)'^, 

R R 

where 

^EiR)A = nA-^-^ctg^SABcn''D^ + ^Da, (134) 

'^E{R)a = nA-^-^ctg^eABcn''D^ -^Da, (135) 

r R 2r 

^E{R)^ = n^dr + Rsin^ -s'^^'^nBdnc + ^ sin— dn^, (136) 

R 2 R 

-E{R)^ = n^dr-Rsin'^^e^'^^nBdnc + ^sin%dn^. (137) 

R 2 R 

The basic Lie brackets are as follows: 

[^E{R)a,^E{R)b] = -^SAB^^E{R)c, (138) 

rE{R)A/E{R)B] = leAB^''EiR)c, (139) 
['EiR)A,''E{R)B] = 0. (140) 

Therefore, the corresponding structure constants are given as follows: 

2 2 
C'(-R) BC = BC, ^C{R) BC = BC, 

and the trivial (cosmetic) Kronecker shift of indices is meant. 
The resulting metric 

g{R) = Sab'E{R)^ (g, 'E{Rf = dAB^'EiR)^ O "-EiRf 
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has the coordinate form 



i.e., 

ds^ = dr^ + sin^ ^ (di?^ + sin^ Mlp^) 
R 

= dr'^ + R^ sm"^ —dn ■ dJT, (141) 
R 

or in more sophisticated terms: 

g{R) = dr®dr + R^ sir? —SABdn^ ® dn^ . 

R 

The algebraic Killing metric ^ab is now given as follows: 

lAB = --^°AB, 



therefore, 



The above analytical formulas may be geometrically interpreted in such a 
way that the unit sphere in R4 is subject to the i?-factor dilatation resulting in 
5*^(0, i?), i.e., the sphere of radius R and origin in R^. The metric g{R) may 
be obtained by restriction of the Kronecker metric in R** to 5'^(0, R). 

Parametrizing S^{Q,R) C R^ by (r, t?, iy9)-coordinates, i.e., 

= R sin — sin d cos Lp^ — R sin — sin -d sin Lp, 

R R 

x^ = Rsin — cosd, x'^ — Rcos — , 

R ' R' 

and expressing 

ds^ = {dx^f + {dx^f + {dx'f + {dx^f 

through (r, t?, iy9), one obtains exactly (|14ip . 

It is seen from ([T34l) . ([T35t . p36)) . p37)) . (fT38l) . p39l) . and (fMO)) that when 
the limit transition i? ^ oo is performed, then ^E{R)a and ^E{R)a become 
asymptotically commutative, then 

lim 'E{R)a = lim '^E{R)a - t^. 

i?— »oo _R^oo OX 

Similarly, 

lim = lim ""EiR)"^ = dx^ , 

lim g{R),j = 5^, 
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i.e., as expected, one obtains the R^-relationships (flat space). 

Now we make use of the pecuharity of the dimension three, i.e., the identifica- 
tion between skew-symmetric second-order tensors and axial vectors. Therefore, 
the matrices of angular velocity will be represented as follows: 



dt 



^diiR)^B = —£^BC^driR)^ , ^^rl'^B — —^TT^'-P B — —S'^BC^rl^ , 



and conversely. 

As always, the capital indices are moved in the "cosmetic" Kronecker sense. 
There is no radius label R at f^ri because these quantities are independent 
of i?. 

After some calculations one obtains the expressions for the basic aholonomic 
velocities ("angular velocities") for the motion on the sphere S^{0,R): 

at at 

where 7c is the rotation vector on S0(3, M) which parameterizes the manifold of 
internal degrees of freedom in the sense of ()131|) with Jc substituted instead of 
k: 

S0(3,M) 3 ip{l<) = cxp (x-^Aj) . 

Denoting the canonical momenta conjugate to r*, respectively by pi, ni, we 
have the expressions for the aholonomic momenta (canonical angular momenta) 
conjugate to Qdr{R}^, ^ri^- 

'Sd,iR)A = 'E\iR;r)p,, ''Sdr{R)A = E' A{R;r)p^, 



Their Poisson brackets are given as follows 

2 
R 



{''Sdr{R)A,'Sdi{R)B} — ■J^^Ab'" ' Sdi{R)c, 



V Sdi{R)A,^ Sdi{R)B} — — t^eas'" ^<5'dr(-R)c 

K 

{'SdAR)A/Sdr{R)B} - 0, 



{'•SriA, '•S'ris} — Sab'" ^Sj 

fr Q r n T ^ C r n 

{ JrlA, '-'rlBl — ~SAB <frlCj 

{^S^\aJS^\b} = 0. 



rlC, 
C r ( 
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The drive and relative quantities are mutually in involution, e.g., 

{'S'dr(-R)A, '•S'ris} = 0, 

and so on. After some calculations we obtain the expression for the kinetic 
energy 

There are a few interesting features of this formula, namely, existing of the mass- 
modifying term J/i?^ in the purely translational part and the characteristic 
gyroscopic term of coupling between two angular velocities. 

The most convenient procedure of deriving equations of motion is one based 
on the Hamiltonian formalism and Poisson brackets. We do almost exclusively 
with the potential models when Lagrangian has the form 

L = T-V{r,K). 

The Legendre transformation expressed in aholonomic terms is as follows: 

_ dL _ dT 

i.e., explicitly, 

^s{r)a = (m + ^yn{R)A-^^nriA, 

The corresponding inverse rule reads 

m mR 

where the raising and lowering of capital indices are meant in the "cosmetic" 
Kronecker-delta sense. 

After substituting all this to the kinetic energy formula, we obtain the ex- 
pression for the kinetic Hamiltonian 

T = -Ls^B Is{R)a'S{R)b + ^6^'' 'S{R)a'Sab 
2m mR 
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When some potential term is admitted, then the Hamiltonian has the form 

H = T + V{r,Ji). 

For the purely geodetic models, when V = Q, the Poisson-bracket form of equa- 
tions of motion, i.e., 

leads to the dynamical system for canonical angular momenta 

^ - ^..-'5„e'5,«,c, (143) 

It is interesting that only the interference term in (|142p contributes to these 
equations because the first and the third terms have vanishing Poisson brackets 
with ^S{R)a: ^S-[\a (they are Casimir invariants). Using the three-dimensional 
vector notation in R^, we can write the above equations as follows: 



d'S{R) 



dt mR^ 
d^S7i 1 



'S,i X 'S{R), (145) 



where there is meant the standard vector product in M.^. 

The Poisson rules quoted above imply that the squared magnitudes of an- 
gular momenta, i.e.. 



S{R)\ = r^'5(i?)A'5(i?)B, 

la |2 xABIq lo 

Jrl — i>i-\A Oris, 



are constants of motion in virtue of the equations above (|143p . (|144p . (|145p . 
p46p . And moreover, these equations imply that the R^-vector 

fs{R)A + 'S,iA (147) 

is a constant of motion. This implies in particular that the R'^-scalar product 
of angular momenta 

<5^^ ^S{R)/S,,B 



and the angle between 'S'(-R), '^ri are constants of motion ftmctionally depen- 
dent on the previous ones. 

In this way the six-dimensional dynamical system (|143p . (|144p . (|145p . (|146p 
has five independent constants of motion and there is only one time-dependent 
degree of freedom in the ^'S'(_R), 'S"^ -space. The only essentially active degree 
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of freedom is the orientation (polar angle) of the two-dimensional plane in 
spanned by the vectors ''S{R), 'Sri and rotating about the axis given by the 

vector ([TT71) . 

Remark: it has been told that the motion of vectors ''S{R), '5^ is indepen- 
dent of the first and second terms in (jl42|) . in particular, of the inertial moment 
/. But of course these terms and the quantity / are relevant for the total de- 
scription of motion. Namely, / appears in the Legendre transformation and its 
inverse. Therefore, the time evolution of configuration variables (r, Jc) depends 
explicitly on / and on all terms in the kinetic energy expression. 
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